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This  technical  note  describes  an  adapted  digital  simulation  program 
for  the  wave  period  processor.  It  was  written  by  Code  D603,  Naval 
Undersea  Warfare  Center,  San  Diego  Division,  Subproject  SF  11  1P1  100. 
Task  11197  (NUWC  Problem  El-19).  This  document  has  been  prepared  because 
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at  the  Naval  Undersea  Warfare  Center  (NUWC)  and  to  a few  persons  outside 
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I.  INTRODUCTION 

^ The  purpose  of  this  technical  note  is  to  document  the 
adaptation  of  a Wave  Period  Processor  (WPP)  Digital  Simulation  Program 
to  NUWC-SD  facilities.  ) 

The  Sperry  Gyroscope  Company  of  Long  Island,  New  York 
developed  a digital  simulation  program  of  the  Wave  Period  Processor 
"to  confirm  the  performance  characteristics  relating  to  the  incentive 
test  measurements".  This  simulation  program  was  made  available  in 
December  1967  to  the  PAIR  Project  Office,  Code  D55^- 

Code  D603,  in  its  quest  to  enhance  in-house  capabilities 
relating  to  digital  simulation  of  signal  processors,  undertook  to 
adapt  and  implement  the  Sperry  Programs  on  the  AN/USQ-20  computers. 

The  adaptation  of  these  programs  to  NUWC  facilities  will  result  in  the 
following  benefits: 

a.  PAIR'S  system  analysis  group  will  be  able  to  analyze  and 
predict  systems  performance  in  a timely  manner. 

b.  In-house  capability  will  be  enhanced  by  exposure  to  the 
approaches  being  used  by  others  who  are  active  in  the  field. 

II.  PROGRAM,  GENERAL 

The  block  diagram  Wave  Period  Processor  digital  simulation  is 
shown  in  Figure  1. 

— The  major  task  in  the  adaptation  of  this  simulation  program 
consisted  of: 

a.  generation  of  a new  random  number  generator. 

b.  .Rewriting  of  the  routines  for  implementation  on  the  AN/uSQ-20 
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The  Wave  Period  Processor  program  is  modularized  into  3 
main  independent  subroutines. 

a.  NOISE  subroutine 

b.  WAVE  subroutine 

c.  PROCS  subroutine 

The  above  subroutines  perform  the  following  functions: 

a.  The  NOISE  Subroutine 

1.  Generates  random  gaus sian-distributed  noise. 

2.  Correlates  the  random  noise. 

3.  Interpolates  and  centers  the  noise  at  20  kHz  for 
the  resulting  narrow  band  gaussian  noise. 

b.  The  WAVE  Subroutine 

1.  Combines  narrow  band  gaussian  noise  and  linear  FM 
signal  according  to  specified  s/N  ratios. 

2.  Measures  each  twelve  period  interval  of  the  wave. 

3.  Computes  the  modulo-32  integer  values  of  the  wave 

measurements. 

c.  The  PROCS  Subroutine 

1.  Generates  the  zone  logic 

2.  Performs  range  bin  detection 

3.  Computes  statistical  information 

In  order  to  better  understand  the  functions  of  the  WPP 
Simulation  Programs  as  related  to  the  hardware,  a brief  description 
of  the  hardware  system  will  be  presented. 
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III.  WAVE  PERIOD  PROCESSOR  DESCRIPTION* 

' 

The  WPP  is  an  active  search  signal  processor.  It  receives 
signals  continuously  from  all  48  active  beams  and  outputs  digital 
amplitude  as  a function  of  bearing  and  time  (range).  The  output 
data  appears  on  the  B-Scan  CRT  display  as  horizontal  strokes  in  a 
matrix  of  48  bearing  bins  by  240  range  bins.  Strokes  may  depict 
returns  from  a single  transmitted  ping  or  from  five  pings,  in  which 
case  the  240  range  bins  are  reduced  to  48,  each  of  which  accommodates 
5 strokes  representing  the  same  range. 

The  48  beams  are  amplified  and  heterodyned  to  20  kHz  at  which 
frequency  very  narrow  bandpass  filters  are  available  (about  440  Hz). 

All  48  channels  are  then  clipped  and  passed  through  "digitizers".  A 
"digitizer"  is  a frequency  measuring  device  which  expresses  the 
frequency  of  the  input  waves  as  an  output  5-bit  number.  With  the  440  Hz 
input  bandwidth  divided  into  32  intervals.,  each  of  the  possible  32 
output  numbers  represents  a frequency  interval  of  some  14  Hz.  Forty- 
eight  zone  logic  and  integrator  circuits  examine  the  modulo-32  numbers 
from  the  digitizers  for  a pattern  of  counts  which  reveals  the  presence 
of  the  linear  FM  slide  in  the  received  signal.  Logic  zones  are 
established,  and  input  signal  counts  are  made  in  these  zones.  Noise  or 
dissimilar  input  signals  produce  equal  counts  in  all  zones;  deviations  are 
indicative  of  signal  returns,  the  magnitudes  representing  the  signal 
strength.  "Integration  and  peak  detection",  as  used  here,  refer  to 
a process  of  accumulating  samples  and  storing  the  largest  counts  for 


* PAIR  System  Docket  9287-02120,  Sperry  Rand  Corporation 
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outputting  serially  through  Display  Control  to  the  B-Scan  CRT.  A count 
is  obtained  on  each  of  the  L8  bearing  channels  every  720  microseconds 
which  is  every  15  microseconds  per  datum. 

IV.  NOISE  SUBROUTINE 

A.  Noise.,  General 

As  the  first  stage  of  the  Wave  Period  Processor  simulation, 
the  noise  simulation  is  the  standard  representation  of  bandpass  noise  as 
the  sum  of  uncorrelated  quadrature  components.  The  random  (pseudo-random) 
noise  generator  provides  two  independent  samples  and  this  uncorrelated 
gaussian  noise  is  passed  into  identical  digital  filters  synthesized 
from  the  impulse  response  of  the  low  pass  equivalent.  U70  Hz  bandpass 
filter  in  the  system.  The  outputs  from  the  filters  are  the  amplitudes 
of  the  in-phase  and  quadrature  components  of  the  correlated  noise. 

These  outputs  are  respectively  multiplied  by  the  sine  and  cosine  of 
the  carrier  frequency  to  obtain  the  in-phase  and  quadrature  components 
of  the  correlated  narrow  band  noise  at  bandpass.  The  components  are 
then  summed  and  recorded  on  magnetic  tape.  Constants  for  the  digital 
filter  used  to  simulate  the  470  Hz  filter  are  generated  by  another 
program. 

B.  Noise,  Detailed 

The  noise  program  KAVEE  and  its  subroutines  BLOCK,  GAUSS, 
and  RANDM  generate  approximately  1 6 seconds  of  simulated  real-time  random 
noise.  Each  of  the  101  blocks  contains  16.128  independent  samples, 
each  sample  or  word  representing  10  microseconds  of  data.  This  block 
size  corresponds  to  the  requirement  of  the  wave  period  program  in  which 
the  signal  wave  repeat  cycle  is  0.16128  seconds. 
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The  random  number  generator  RANDM  produces  a pseudo-random 
number  uniformly  distributed  between  zero  and  one.  The  routine  GAUSS'* 
obtains  two  random  numbers  from  RANDM.,  and  R^  and  in  turn  delivers 
G1  ■ (-2  In  3in  2rTRg  and  Gg  = (“2  In  R^)1^2  cos  2rrR2  which  are 

independently  gaussian  distributed  with  zero  mean  and  unit  variance. 

The  subroutine  BLOCK  performs  the  major  computations  for  one 
block  of  noise  (l6l28  samples).  Input  parameters  to  the  noise  program 
for  sampling  time  (t).,  Bandwidth  (BW),  and  the  Z-transform  coefficients 
(C1-C9)  of  the  Butterworth  filter  are  computed  elsewhere  by  the  program 
CONSTS  (plus  routines  BUTRW5  and  INVRT6)  and  are  read  from  cards  by 
BLOCK.  These  are  read  only  once.,  since  the  same  values  are  used  for 
each  of  the  101  blocks  of  noise  which  comprise  one  complete  noise  run. 

For  the  first  runs,  T=50  microseconds.,  BW=1+71.2  Hz. 

The  two  gaussian  random  numbers  represent  uncorrelated 
noise  samples  for  the  sine  and  cosine  channels  with  zero  mean  and  unit 
variance.  These  are  produced  every  50  microseconds.  The  uncorrelated 
samples  are  passed  through  the  two-step  recurrence  equations  to  obtain 
low-pass  correlated  noise. 

For  Sine  Channel: 

YSi  = c6  • xs1_1  + c7  • xsi_2  + Cg  • xsi-5  + C9  • xst_k 

where : • 

XSi  = C1  • xsi_1  - c2xsi_2  + c5  • xsi-5  - C4  • xsi-4  + c5  • xsi_5  + 


Complete  descriptions  of  these  routines  are  enclosed  in  Appendix  A 
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For  Cosine  Channel; 

IC1  ' C6  ' xol-l  * °T  ' x0l-2  ♦ °8  ’ XCi-3  * °9  • xcl-l. 
where; 

XC1  • C1  • XC1-1  • °2  • xcl-2  * °3  ■ xol-3  - ck  • xcl-4  * C5  • xci-5  * °2 
C1-C9  are  the  Z-transfora  coefficients  of  the  Butterworth  filter.  G ^ is 
the  zero  mean,  unit  variance  random  number  generated  by  subroutine  GAUSS. 

A three-point  Lagrange  interpolation  is  applied  to  the  low- pass 

outputs  of  the  second  stage  recurrence  equations  to  produce  one  sample 

every  10  microseconds.  Narrow  band  noise  is  desired,  centered  at  2C  kHz. 

The  low-pass  outputs  are  therefore  multiplied  respectively  by  the  sines 

and  cosines  of  multiples  of  72°,  which  correspond  to  10  microseconds  of 

a 20  kHz  wave,  and  the  two  channels  are  combined. 

For  wl  = sin  lt)  t.  and  w2  = cos  'u  t.  . 

o l o 1 1 

= Lx  * wl  • YSi_2  + L2  • wl  • YSi_1  + L?  • wl  • YSi  + Lj^  * w2  • YC^ 

+ Lp  • w2  • YCi_1  + L,  • w2  . YC^ 

where  is  the  Lagrange  interpolation  constant,  YS^'s  and  YC^'s  are 
the  low-pass  outputs  from  the  sine  and  cosine  channel  second  stage 
recurrence  equations,  is  the  desired  channel  band  center,  and  Y^  is 
the  final  output  value. 

In  actuality,  the  constants  L^wl  through  L^w2  in  the  above 
equation  are  pre-multiplied.  The  L^'s  are  the  three-point  interpolation 
constants  for  the  intervals  -0.4,  -0.2,  0.0,  +0.2,  +0.4.  The  w^'s  are 
the  numerical  values  for  the  sines  and  cosines  of  72°,  36°,  and  0°. 


p.70,  Lagrangian  Methods,  Introduction  to  Numerical  Analysis,  F.B.  Hildebrand 
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Five  noise  samples,  Y^,  are  therefore  generated  each  cycle  of 
the  recurrence  and  interpolation/centering  equations,  one  each  10 
microseconds  for  the  original  50  microsecond  sample.  However,  since  5 
is  an  uneven  factor  of  16128,  the  noise  block  size,  a modulo-  5 indexing 
scheme  has  been  employed.  In  a set  of  5 blocks: 

J = 0:  Block  1 - l6lJ0  samples  are  computed. 

J = 2:  Block  2-2  extra  samples  from  first  block  begin  second 
block,  and  16130  new  samples  are  computed. 

J = 4;  Block  3-4  extra  samples  from  second  block  begin  third 
block,  and  l6l25  new  samples  are  computed. 

J = 1:  Block  4-1  extra  sample  from  third  block  begins  fourth 
block,  and  16130  new  samples  are  computed. 

J = 3:  Block  5-3  extra  samples  from  fourth  block  begin  fifth 
block,  and  l6l25  new  samples  are  computed. 

The  extra  samples  computed  for  one  block  of  noise  are  moved 
to  the  beginning  of  the  next  block  and  the  noise  function  is  therefore 
continuous.  A maximum  of  16132  samples  may  be  present  at  one  time,  even 
though  only  16128  samples  are  used  in  any  one  noise  block  for  output. 

The  subroutine  BLOCK  concludes  by  writing  one  block  of  16128 
narrow  band  gaussian  noise  values  on  magnetic  tape  and  computing  the 
mean  of  the  sum  and  the  mean  of  the  sum  of  the  squares  for  the  l6l28 
noise  values. 

The  control  routine  for  the  noise  generation,  KAVEE,  calls 
on  the  subroutine  BLOCK  a total  of  101  times  for  one  complete  noise  run. 
Each  time  the  control  routine  calls  the  subroutine  BLOCK,  plus  once  after 
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the  final  call,  it  will  print  the  BLOCK  input  parameters,  including  the 
current  block  number,  indices,  modulo-  5 index,  initial  entry  flag, 
current  starting  point  of  the  random  number  generator,  and  current 
values  of  the  two  step  recurrence  equations.  The  control  routine  will 
output  the  mean  and  mean  of  squares  values  and  compute  and  output  the 
varian-e  and  standard  deviation,  all  for  each  noise  block.  At  the  end 
of  a complete  run,  it  will  also  compute  and  output  the  mean,  mean  of 
squares,  variance,  and  standard  deviation  from  the  mean  values  of  all 
noise  blocks  combined. 

The  total  output  for  101  blocks  of  noise  is  approximately 
1500  feet  of  magnetic  tape.  Since  each  group  of  16,128  words  is  needed 
to  generate  22k  samples  of  the  wave  period  measurements,  one  complete 
noise  tape  will  provide  101  x 22k  = 22,62k  samples  of  wave  period 
measurement. 

Symbol  Correspondence 

ABS  Library  routine  for  absolute  value  of  a real  argument. 

ALOG  Library  routine  for  natural  logarithm  of  a real  argument. 

AHOLE  Second  of  the  pair  of  random  numbers  used  to  compute  two 

gaussian  distributed  numbers. 

BLKS  Maximum  number  of  noise  blocks,  in  floating  point.* 

BLOCK  Name  of  subroutine  which  performs  the  noise  computations. 

BW  Bandwidth  of  Butterworth  filter  (L7I.2  Ho).** 


*Mixed  mode  arithmetic  is  not  available;  therefore,  constants  must  be 
entered  in  any  and  all  modes  in  which  they  are  used. 

•♦Numbers  in  parentheses  refer  to  values  used  in  first  NUWC  run. 
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C1-C9 

CA1-CA5 

CB1-CB3 

COS 

DENOM 

DUMMY 

G 

GAUSS 

II 

IF  LAG 


IMAX 

I RAJ©  3 

J 


♦Refers 


Z-transform  coefficients  of  Butterworth  filter. 

Three  point  Lagrange  interpolation  constants  for  the 
cosine  channel. 

Library  routine  for  trigonometric  cosine  of  a real  argument. 
Number  of  samples  for  one  noise  block.,  in  floating  point. 

Dummy  argument  used  to  set  the  random  number  generator. 

Array  containing  l/MULTP  of  a noise  block  (&88). 

(previously*  equal  to  length  of  one  magnetic  tape  record) 
Function  subroutine  which  computes  gaussian  distributed 
random  numbers. 

Indicator  for  number  of  noise  blocks  written  on  output  tape. 
Indicator  for  first  or  second  entry  to  routine  GAUSS 
IFLAG  = 0.  pick  up  sine  channel  gaussian  value 
IFLAG  = 1,  pick  up  cosine  channel  gaussian  'value  (no  longer 
used  in  routine  BLOCK) 

Major  processing  index  for  the  5_step  noise  generation 
equations. 

Entry  to  random  number  generator  that  will  obtain  the  last 
random  number.,  which  is  also  th?  starting  point  for  the  next 
random  number. 

Modulo  5 counter  for  shifting  extra  values  (>  16128)  from 
the  end  of  one  noise  block  to  the  beginning  of  the  next  block. 


to  Sperry  program  format. 
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KAVEE 

LNGTH 


MLTM1 

MOD 

MULTP 


NBLKS 

NBLOCK 

NTRY 


NUMBR 


RANDM1 


RANDM2 


RNDSET 


Main  control  routine  for  noise  generation. 

Intermediate  number  of  noise  points  generated  equalling 
l/MULTP  of  cne  block  of  noise  (288)  (originally  equal  to 
number  of  noise  samples  on  one  magnetic  tape  record.*) 

Equal  to  MULTP-1  for  indexing  purposes  (55). 

Library  routine  for  remaindering  integer  arguments. 
Multiplication  factor  which,  along  with  LNGTH,  determines  the 
length  of  1 block  of  noise  (56)  (previously  equal  to 
the  number  of  magnetic  tape  records  for  one  noise  block.*) 
Maximum  number  of  noise  blocks. 

Index  and  indicator  for  current  noise  block  being  processed. 
Flag  which  indicates  initial  pass  in  a given  execution  of 
program. 

= 0,  initial  pass 
= 1,  succeeding  pass 

Number  of  noise  samples  in  one  noise  block 
= MULTP  x LNGTH 

Entry  to  random  number  generator  to  obtain  the  next  pseudo- 
random number. 

Entry  to  random  number  generator  to  store  a new  starting 
point  for  generator. 

Octal  number  used  to  start  the  random  number  generator  at 
a particular  point. 


♦Refers  to  Sperry  program  format. 
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MLTMl  Equal  to  MULTP-1  for  indexing  purposes  (55) 

MOD  Librarv  routine  for  x.  _ 
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SA1  -SA.' 

SB1-SB' 

SIOY 

SIG2 

SIN 

SQRT 

SS4Y 

SSQ”. 


SUMY 

SUMZ 

T 

TKMP 

THANK 

VARY 

var:i 

xc 

xs 


Y 

YBAK 

YC 


Three  point  Lagrange  Interpolation  constants  for  the  sine 
channel 

Standard  deviation  for  one  block  of  noise. 

Standard  deviation  for  all  blocks  of  noise  combined. 
Library  routine  for  trigonometric  sine  of  a real  argument. 
Library  routine  for  square  root  of  a real  argument. 

Sum  of  the  squares  of  samples  for  one  block  of  noise. 

Sum  of  the  squares  of  samples  for  all  blocks  of  noise 
combined. 

Rum  of  samples  for  one  block  of  noise. 

Ruin  of  samples  for  all  blocks  of  noise  combined. 
Butteivorth  filter  sampling  time. 

Temporary  location  for  Intermediate  calculations . 

First  of  the  pair  of  random  numbers  used  to  compute  two 
gauss  inn  distributed  numbers. 

Variance  of  samples  for  one  block  of  noise. 

Variance  of  samples  for  all  blocks  of  noise  combined. 

Array  containing  the  previous  four  and  the  current  first- 
step  recurrence  noise  numbers  for  the  cosine  channel. 

Array  containing  t he  previous  four  and  the  current  first- 
step  recurrence  noise  numbers  for  the  sine  channel. 

Array  containing  noise  samples. 

Mean  of  samples  for  one  block  of  noise. 

Array  containing  the  previous  two  and  the  current  second- 
step  recurrence  noise  numbers  for  the  cosine  channel. 
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YMSQ 

Mean  ■ 

YS 

Array 

step 

ZBAR 

Mean 

ZMSQ 

Mean 
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INPUT 


U punched  cards  contain  input  values  for  the  program: 


Card 

Variable 

Columns 

Format* 

1 

T 

5-27 

E23-16 

BW 

32-54 

E23.16 

2 

C1 

4-26 

E23.16 

C2 

31-53 

E23.16 

C3 

58-80 

E23.16 

3 

c4 

5-27 

E23.16 

C5 

32-54 

E23-16 

4 

c6 

5-19 

E15.8 

cT 

24-3  8 

E15.8 

c8 

43-57 

E15.8 

CQ 

62-76 

E15-8 

In  the  case  of  the  noise  program,  double  precision  variables  T.  BW., 

C1-C5  are  entered  through  an  E field. 

OUTPUT 

A list  of  the  ^.-transform  coefficients,  the  sampling  time,  and 
bandwidth  of  the  Butterworth  filter  are  printed  out  as  soon  as  they  are 
read  in  from  the  cards. 

At  the  beginning  of  each  noise  block  computation,  the  noise 
block  number,  two  indices  which  combined  indicate  the  length  of  a noise 
block,  the  modular  5 index,  initial  pass  flag,  and  the  last  number  (in 
octal  integer  form)  generated  by  the  random  number  routine  are  printed. 


#"In  format  statements  used  for  data  input,  F,.,  F,.  and  G field  specifications 
are  interchangeable.  In  addition,  Double  Precision  variables  may  be 
entered  with  those  specifications".  Fortran  IV  for  the  NUWO  (SP)  1130 
and  AN/USQ-20  computers. 
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The  last  random  number  Is  useful  if  one  wishes  to  stop  the  program 
and  restart  it  later  at  the  point  the  random  number  generator  left 
off.  The  call,  RANDM2  (RNDSET),  to  the  random  number  generator  will 
put  the  variable  RNDSET  in  as  the  starting  point. 

The  sine  and  cosine  recurrence  tables  XS,  XC,  YS,  YC  are 
printed  out  so  that  they  may  also  be  entered  to  restart  the  program.* 

At  the  end  of  each  noise  block  computation,  the  mean  of 
the  sum,  mean  of  the  sum  of  squares,  variance,  and  standard  deviation 
of  the  l6l28  samples  are  printed  out. 

After  one  entire  noise  run  (101  blocks)  has  been  completed, 
the  mean  of  the  sum,  mean  of  the  sum  of  squares,  variance,  and 
standard  deviation  are  printed  for  all  noise  blocks  combined. 

The  101  noise  blocks  of  l6l28  samples  each  are  written 
identically  on  two  IBM  compatible  magnetic  tapes,  and  an  end  of  file 
mark  completes  the  output  on  each  tape. 

Subroutines 

Descriptions  of  functions  and  entry  points  of  the  routines 
GAUSS  and  RANDM  follow. 

Library  functions  used  include: 

MOD  Computes  remainder  of  two  integer  arguments. 

SQRT  Computes  square  root  of  a real  argument, 
plus  those  called  by  the  routine  GAUSS 

SQRT  (Above) 

ABS  Computes  absolute  value  of  a real  argument. 


* See  "Consecutive  Runs  of  Noise  Program" 
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ALOG  Computes  natural  logarithm  of  a real  argument. 

SIN  Computes  trigonometric  sine  of  a real  argument. 

COS  Computes  trigonometric  cosine  of  a real  argument. 

Execution 

Fortran  IV  on  the  NUWC  (SD)  12.50  and  AN/USQ-20  computers 
is  an  automatic  compile  and  execute  system,  so  the  entire  program  deck, 
as  listed  following,  is  used  for  a run.  A PAUSE  has  been  inserted 
to  stop  the  program  at  the  beginning  of  execution  for  the  purpose  of 
mounting  the  output  tapes. 

Tape  unit  assignments  for  NUWC-SD  Univac  1250  or  AN/USQ-20 

computers: 

System  tape  on  Unit  Ml 
Noise  output  tape  on  Unit  M2 
Noise  output  tape  on  Unit  M3 
Execution  Time 

To  generate  one  block  of  noise  data  of  l6l?8  samples  takes 
25  minutes  on  the  Univac  AN/USQ-20  computer  and  1 minute  15  seconds 
on  the  Univac  1230  computer.  However,  in  order  to  take  advantage  of 
the  fast  double  precision  hardware  on  the  1230  computer,  the  entire 
program  was  converted  to  double  precision,  with  a running  time  of  50 
seconds  per  noise  block. 

Accuracy 

Double  precision  arithmetic  is  used  in  the  recurrence 
formulas  to  generate  the  XS^  and  XC^  values.  Single-precision 
arithmetic  is  used  to  generate  the  noise  values  for  output  on  magnetic 
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tape.  However,  as  stated  above,,  the  entire  program  was  also  written 


in  double  precision  to  take  advantage  of  the  fast  floating  point 


hardware  on  the  Univac  1230  computer.  The  program  in  both  forms  is 


included  in  Appendix  C 


Consecutive  Runs  of  Noise  Program 


When  more  than  101  blocks  of  continuous  noise  are  desired 


the  last  values  from  the  two-step  recurrence  equation  tables  and  the 


last  starting  point  of  the  random  number  generator  must  be  initially- 


entered  into  the  noise  program.  Seven  cards  contain  the  generator 


starting  point,,  RNDSET.  the  five  values  for  the  first-step  recurrence 
equation  for  both  the  sine  (XS)  and  cosine  (XC)  channels,  and  the  three 
values  for  the  second-step  recurrence  equation  for  the  sine  (YS)  and 
cosine  (YC)  channels.  The  double  precision  noise  pregram  included  in 


this  document  contains  the  instructions  for  loading  the  above  values 


CARD 


VARIABLE 


COLUMNS 


FORMAT 


1 RNDSET  9-20  012 

2 XS(l)  4 - 26  D23.16 

XS(2)  30  - 52  D23.16 

XS(3)  56  - 78  D23.16 

3 XS(4)  4-26  D23.16 

XS(5)  50  - 52  D23.16 

4 XC(1)  4-26  D23.16 

XC(2)  30  - 52  D23.16 

XC(3)  56  - 78  D23.16 

5 XC(4)  : 4-26  D23.16 

XC(5)  50  - 52  D23.16 

6 YS(l)  4-26  D23.8* 

YS(2)  30  - 52  D23.8 

YS(3)  56  - 78  D23.8 

7 YC(l)  4-26  D23.8 

YC(2)  30  - 52  D23.8 

YC(3)  56  - 78  D23.8 

*The  second-stage  values  would  be  entered  through  E w.d  field  for 
the  primarily  single  precision  noise  program. 
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WAVE  SUBROUTINE 


A.  Wave,  General 

The  WAVE  program  is  the  second  stage  of  the  simulation, 
using  as  input  the  output,  on  magnetic  tape,  from  the  NOISE  program. 

The  noise  samples  are  combined  with  a linear  FM  slide  which  is  scaled 
to  provide  the  proper  input  signal-to-noise  ratio.  The  combined 
signal  plus  noise  samples  are  then  observed  for  a sign  change  from 
minus  to  plus.  To  increase  the  accuracy  of  the  wave  period  measurement, 
an  interpolation  is  made  for  the  first  and  thirteenth  zero  crossing 
and  for  every  thirteenth  zero  crossing  thereafter  to  determine  the 
exact  time  of  zero  crossover.  The  time  for  thirteen  of  these 
crossings  denotes  twelve  periods  of  the  wave  and  this  time  is 
presented  to  a mod-32  counter  driven  by  a simulated  digital  clock  at 
a 2.424  megacycle  rate.  The  output  of  this  counter  is  a number 
between  0 and  31  and  occurs  every  720  microseconds  of  the  data, 
producing  224  wave  measurements  every  0.16128  seconds.  The  integer 
mod-32  numbers  are  written  on  magnetic  tape  for  input  to  the  third 
stage  of  the  WPP. 

B.  Wave,  Detailed 

The  wave  period  measurement  program,  WAVE,  generates  a 
signal  wave  with  a repeat  cycle  of  0.16128  seconds.  Gaussian  random 
noise  is  combined  with  the  signal  according  to  specified  input  signal- 
to-noise  ratios.  A twelve  period  interval  of  the  resulting  wave  is 
measured  every  720  microseconds  , truncated  to  integers  , taken  mod-32 , 
and  written  as  output.  This  produces  224  measurements  per  cycle  for 
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each  S/N  ratio,  of  which  12  were  used  in  the  first  run  of  the  WAVE 
program.  The  result  is  224  x 12  = 2688  measurements  per  output 
record  and,  considering  101  cycles,  or  101  distinct  noise  blocks, 

101  x 2688  = 271,488  total  wave  measurements  for  one  run  of  the  WAVE 
program. 

The  16  seconds  of  narrow  band  gaussian  noise  from  the  NOISE 
program  provides  the  input  to  the  wave  period  measurement  simulation 
program.  The  WAVE  program  adds  this  noise  data  to  different  levels 
of  signals  to  generate  the  signal  plus  noise  waves  at  signal- to-noise 
ratios  of  -®,  -6,  -4,  -3,  -2,  -1,  0,  +1,  +2,  +3,  +4,  and  +6  db  for 
the  first  run  of  WAVE,  the  first  input  s/N  ratio  being  noise  only. 

The  signal  is  a linear  FM  slide  generated  as  follows: 

A frequency  modulated  signal  can  be  expressed 

a(t)  = A cos  0 (t), 

where  the  amplitude  A is  kept  constant  and  the  instantaneous  frequency 
of  the  cosine  function  is  varied  in  accordance  to  the  signal. 
Instantaneous  frequency  is  defined  as 

2tt  ., 

which,  for  a linear  FM  slide,  may  be  expressed  as 

F + AFt 


where 


F = carrier  frequency  ='  20,000  hz 
AF  = peak  frequency  deviation  per  slide  time 
Integrating  the  above  expression  with  respect  to  time  yields 
0(t)  = 2tt  J(F  + AFt)dt  = 2n(Ft  + AFt£>/2)+eo 
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For  the  WPP,  the  peak  frequency  deviation  is  computed  as 

. _ Frequency  bandwidth 440  cycles/sec,  _ q , / 

Slide  transmission  time  - 0.16128  sec.  ~ 

Therefore,  assuming  unity  amplitude  and  zero  phase,  the  linear 

FM  signal  can  be  represented  as 

(1)  S = cos  2rr{Ft  + AFt2/2)  = cos  2rr(20,000t  + 1364t2) 

The  signal  has  a repeat  cycle  of  0.16128  seconds  or  161,280 

microseconds.  A sample  time  of  10  microseconds  was  chosen.  This 

means  that  the  strength  of  the  signal  (and  the  random  noise)  at 

the  end  of  every  10  microsecond  interval,  i.e. , at  16128  points, 

must  be  calculated.  As  input  to  the  program,  there  are  101  blocks 

of  16128  points  of  white  noise,  allowing  the  generation  of  101 

different  signal  plus  noise  waves  from  the  one  0.16128  second  signal 

wave.  This  is  done  for  each  of  the  12  signal  to  noise  ratios 

specified. 

In  the  generation  of  the  signal  with  equation  (l)  above,  each  signal 
point  J may  be  computed 

S(j)  = cos  2n(20,000t  + 1364t2) 

where 


t = 10  Usee  (j  - 


16129 


) 


Initially  in  the  WAVE  program,  a set  of  constants  are 
computed  for  amplitudes  of  1 db,  2 db , 3 db , 4 db,  6 db,  and  others 
not  used  in  the  first  WAVE  run.  A set  of  12  scale  factors  L are 
computed  for  the  12  S/N  ratios  desired,  using  the  db  constants,  where 
SCALE(L)  = RMSx,/~2  x DB  Constant 
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RMS  is  the  root  mean  square  value  of  the  noise  computed  during 
generation  of  the  noise  and  is  injected  as  a correction  factor. 

The  noise  and  the  signal  are  then  combined  in  the  following  way: 

SPN(J)  = NOISE(J)  + SCALE(L)  x SIGNAL(j) 
for  signal-to-noise  ratio  L. 

The  signal  may  be  combined  with  the  noise  only  at  specified 
Intervals  to  obtain  a more  realistic  waveform.  An  explanation  may 
be  given  by  example  of  the  first  program  run  at  NUWC.  Signal  was 
combined  with  noise  for  each  third  noise  block  so  that,  for  each 
S/N  ratio,  two  noise  blocks  (16128  samples  per  block ) were 
interspersed  between  each  signal  plus  noise  block.  Therefore,  for 
each  noise  or  signal  plus  noise  block 

SPN 

where 


I = S/N  ratio  from  1 through  12, 

remembering  S/N  ratio  (l)  = -®  or  noise  only,  the  following  arrangement 
is  generated: 


S/N  Ratic 

* 


1 

2 

5 

4 

PN1 

spp 

SPN! 

spp 

PN2 

spp 

SPNi 

SPN, 

pp 

SPN1 

spp 

SPN, 

- 99  100  pi 


SFp 

SPN1 

SPN 

SPNp 

The  actual  wave  period  measurements  are  made  by  examining 
the  signal  plus  noise  samples  of  each  above  wave  separately  for  a 
change  in  sign  from  minus  to  plus.  The  simulation  detects  13 
successive  positive  zero  crossings  and,  if  the  first  and  thirteenth 
crossing  points  are  not  identically  zero.,  an  interpolation  is  made 
on  those  two  points,  assuming  a sine  wave, to  obtain  more  accurate 
crossing  times.  The  time  elapsed  between  the  first  and  thirteenth 
crossings  determines  12  periods  of  the  signal  plus  noise  wave.  720 
microseconds  (or  72  - 10  microsecond  samples)  are  added  to  the  first 
crossing  time  and  the  above  process  is  repeated.  Thus,  224  elapsed 
time  measurements  (called  TT's  in  the  program)  will  be  obtained 
every  0.16128  seconds. 

Each  of  the  elapsed  time  measurements  are  then  subjected 
to  a modular  32  counter  driven  by  a simulated  2.424  megacycle  clock. 
Also  injected  is  a sliding  preset  which  is  set  to  zero  at  the  start 
of  each  224  FM  slide  transmission  and  is  preset  to  the  next  higher 
counter  state  every  seventh  sample  period  (i.e.,  every  7 TT's). 

In  an  actual  system,  however,  the  counter  may  start  in  any  arbitrary 
position.  This  effect  is  remedied  by  an  averaging  process  in  the 
third  program.  The  elapsed  time  measurements  (TT’s)  are  therefore 
combined  with  the  clock  and  the  sliding  preset  in  the  manner 

TT  x CLOCK  + PRESET., 

truncated  to  integers.,  and  taken  modulo-32.  The  resulting  integer 
numbers  between  0 and  31  (denoted  IT's  by  the  program  and  having  a 
one  to  one  correspondence  with  the  TT's)  are  written  as  output  on 
magnetic  tape. 
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ibol  Corresuondence 


Library  routine  for  trigonometric  cosine  of  a real  argument 


used  in  various  calculations:  also 


CELT 


10  microsecond  constant 


DBIPHF  Constant  value  for  1 l/2  db  (not  used  in  first  run) 
FO  Carrier  frequency  (20,000  kHz)* 


One-half  of  the  peak  frequency  deviation: 
l/2(bandv id th/cycle  time)  = l/2(U40.0/0. 16128) 


Constant  value  for  3 db 


Constant  value  for  k db 


Constant  value  for  5 db  (not  used  in  first  run) 
Constant  value  for  6 db 


Constant  value  for  12  db  (not  used  in  first  run) 
Constant  value  for  18  db  (not  used  in  first  run) 
Constant  value  for  2U  db  (not  used  in  first  run) 
Constant  value  for  30  db  (not  used  in  first  run) 
Signal  bandwidth  (W+0  Hz) 


F18DB 


F24DB 


F30DB 


Floating  point  value  **  for  the  number  of  positive  zero 
crossings  for  one  wave  measurement  (number  of  wave  periods 
for  each  measurement) 


FNUMBR  Floating  point  value  of  number  of  noise  samples  in  one 


block;  also  the  slide  transmission  time  X 10  or  the  repeat 


* Numbers  in  parenthesis  refer  to  values  used  in  first  NUWC  run. 

**  Mixed  mode  arithmetic  is  not  available;  therefore,  constants  must 
be  entered  in  any  and  all  modes  in  which  they  are  used. 
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ICROSS 

INI’ 


IT 


ISET 


J 


L2,L6 


LNGTH 


LOUT 

LSCALE 


lti 


Index  for  NP  positive  zero  crossings  (l?) 

Library  routine  for  truncation  of  a real  argument  to  an 
integer 

Modulo  - 32  integer  wave  period  measurements  computed  from 
TT's 

Sliding  preset  combined  with  TT's  in  forming  modulo-  32 
IT's 

General  index;  also.,  current  position  within  noise  block 
when  computing  positive  zero  crossings 
Previously  * were  status  locations  for  NTRAN  operations 
which  wrote  IT's  on  magnetic  tape;  Equal  to  LTT  for 
successful  operation.  No  longer  used. 

Intermediate  number  of  noise  points  generated  equalling 
- ^-7 ■ of  one  block  of  noise  (288).  Originally  * equal 
to  number  of  noise  samples  on  one  magnetic  tape  record. 
Number  of  wave  period  measurements  per  cycle  (224). 

Number  of  signal/noise  ratios  to  consider;  therefore, 
number  of  scale  factors  to  calculate  and  the  number  of  S/N 
ratios  to  calculate  for  each  noise  block  (12). 

Number  of  wave  measurements  per  cycle  X number  of  S/N 
ratios  (224  X 12  = 2688);  Size  of  one  magnetic  tape 
record  of  IT's 


* Refers  to  Sperry  program  format. 
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MLTM1  Equal  to  MULTP-1  for  indexing  purposes  (55) 

MOD  Library  routine  for  remaindering  integer  arguments 

MULTP  Multiplication  factor  which,  along  with  LNGTH,  determines 
the  length  of  one  block  of  noise  (56) 

NBLKS  Maximum  number  of  noise  blocks  to  be  used  for  one  WAVE  run 

NP  Number  of  positive  zero  crossings  for  one  wave  period 

measurement  (12) 

NUMBR  Number  of  samples  in  one  noise  block  (l6l28) 

ONEDB  Constant  value  for  1 db 

OSCL  Digital  clock  frequency  calculation 

g 

(2.424  x 10  cycles)  = (2.424  megacycles) 

PI  n( 3. 1415927) 

PLIN  Intermediate  calculation  for  interpolation  of  positive 

zero  crossing;  Equal  to  0+value/(0+ value  - 0_ value) 

RMS  Root  mean  square  of  noise  used  as  a correction  factor 

for  s/N  ratio  scale  factors;  Equal  to  the  standard 
deviation  for  the  total  noise  run  (see  SIGZ  in  program 
KAVEE) 

S Signal  value 

SCALE  Array  containing  scale  factors  for  Signal  + Noise  computations 

SJ  Floating  point  value  for  J when  used  as  the  position  index 

within  the  noise  block 

SPN  Temporary  storage  for  one  block  of  noise;  also  signal  is 

combined  with  the  noise  for  signal  + noise  values  here 
SQRT2  Numerical  value  for  the  square  root  of  2.0  (l. 4142136) 

START  First  positive  zero  crossing  and  therefore  the  start  of 

one  set  of  224  wave  period  measurements 
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Input 

One  IBM  compatible  magnetic  tape  with  101  blocks  of  low 
bandpass  random  noise,  16128  samples  in  each  block. 

Output 

Each  noise  block  input  is  processed  with  each  of  the  12 
S/N  ratios,  and  one  output  wave  record  consists  of  224  x 12  = 2688 
mod-52  integer  wave  period  measurements.  There  will  be  101  such 
2688  word  records  written  on  an  IBM  compatible  magnetic  tape. 

After  each  record  is  written,  the  index  (l-lOl)  is  printed 

out. 

Subroutines 

Library  functions  called  are: 

COS  Computes  trigonometric  cosine  of  a real  argument 
INT  Truncates  a real  argument  to  an  integer 
MOD  Computes  remainder  of  two  integer  arguments 

Execution 

The  entire  program  deck  is  loaded,  since  Fortran  IV  at 
NUWC-SD  is  an  automatic  compile  and  execute  system.  A listing  of 
the  WAVE  program  is  enclosed.  A PAUSE  has  been  inserted  to  stop 
the  program  at  the  beginning  of  execution  so  that  the  noise  input 
tape  and  the  output  tape  may  be  mounted. 

Tape  Assignments  for  NUWC-SD  Uni vac  1250  or  AN/uSQ-20  computers 
System  tape  on  Unit  Ml 
Noise  input  tape  on  Unit  M2 
Scratch  tape  on  Unit  M3 
Wave  output  tape  on  Unit  M4 


CQWWS&pti  i 
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Execution  Time 

To  generate  one  record  of  wave  data  of  2688  words  takes  2 l/2 
hours  on  the  Univac  AN/USQ-20  computer  and  5 minutes  on  the  Univac 
1230  computer.  However.,  in  order  to  take  advantage  of  the  fast 
double  precision  hardware  on  the  1230  computer,  the  entire  program 
was  converted  to  double  precision,  with  an  average  running  time  of 
3 minutes  per  record. 

Accuracy 


The  wave  program  was  written  in  single  precision  since 
the  final  output  is  integer  and  requires  no  greater  precision. 
However.,  as  stated  above.,  the  entire  program  was  also  written  in 
double  precision  to  take  advantage  of  the  fast  floating  point 
hardware  on  the  Univac  1230  computer.  The  program  in  both  forms  is 
included  in  Appendix  C. 

VI.  PROCS  SUBROUTINE 

A.  PROCS.,  General 

The  PROCS  Routine  commences,  as  related  to  hardware, 
with  the  zone  logic  circuits.  See  Figure  2.  The  function  of  the 
zone  logic  circuits  is  to  optimize  the  processing.  In  actual 
practice  the  zone  logic  circuits  function  as  band  pass  filters. 
Each  zone  is  200  Hz  wide  and  they  are  displaced  by  100  Hz. 


f 
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Figure  2 . ZONE  LOGIC  SPECTRUM 
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The  input  to  the  zone  gate  logic  block  are  numbers  0 to  31.  These 
numbers  represent  a frequency  and  in  the  program  are  denoted  as  IT's. 
The  zone  logic  circuits  are  simulated  by  dividing  the  numbers  0-31  into 
zone  A.  B.,  A and  B.  An  "A"  channel  is  generated  for  an  IT  between 
0-15  and  a B channel  for  an  IT  between  8-23  Inclusive.  A moving  sum 
of  224  samples  is  taken  for  each  channel  and  Its  complement,  l.e. , 

sum^  = •£  ZONCT 


sum^  = sumi_1  + ZNC0UNT1+225  - ZONCT^ 


This  moving  sum  is  denoted  as  AHUM  and  its  length  (224  samples) 
represents  the  signal  period. 

As  stated  on  the  description  of  the  hardware  system,  the 
digitizer  is  preset  in  such  a manner  that  a coherent  signal  will 
appear  as  stationary  frequency.  Therefore.,  ANUM  represents  the 
number  of  times  a particular  frequency  falls  on  a zone  and  is  a 
measure  of  period  correlation  during  a signal  period.  The  ANUMS 
for  channels  A and  B are  compared  and  the  maximum  selected  to  denote 
the  ANUM  of  that  particular  sample.  The  ANUM's  are  then  range 
peaked  detected  by  selecting  the  maximum  ANUM  every  6?2  (range  bin) 
ANUMS  in  a specified  range  bin.  In  the  hardware  system  these  ANUMS 
are  the  output  of  the  Wl’P.  In  the  simulation  program  these  numbers 
are  used  to  plot  probability  of  exceeding  threshold  vs  threshold 
with  input  S/N  as  a parameter. 

B.  PROCS,  Petal] ed 

The  PROCS  subroutine  is  divided  into  four  phases. 
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Phase  1 

The  PROCS  routine  begins  by  reading  the  number  of  NRUNS, 
LSCALE,  and  IPEAK;  then  it  positions  the  tape  to  the  file  specified 
by  the  LSCALE.  The  input  to  this  routine  are  numbers  called  IT's 
and  are  stored  by  files  on  magnetic  tape.  Each  file  represents  a 
particular  input  signal  to  noise  ratio.  Due  to  the  core  limitation 
only  11,424  of  the  IT's  are  read  at  first.  The  subroutine  ZONEIT 
is  called  and  each  IT  is  ZONED. 

The  subroutine  ZONEIT  defines  8 separate  zoning  conditions 
where  the  zones  are  shifted  to  cover  one  more  number  to  the  right  and 
one  less  number  to  the  left,  and  where  zones  A,  A,  B,  B each  time  correspond 
to  a different  set  of  numbers  between  0-31.  A sum  of  the  number  of  times 
it  falls  into  a zone  is  kept.  These  sums.,  one  for  Zone  A and  one  for 
Zone  B are  divided  by  8 and  the  results  denoted  as  AIT  or  BAIT. 

The  output  for  the  ZONEIT  subroutine  is  denoted  AIT  and 
BAIT  for  each  IT,  and  these  are  recorded  in  records  of  224  on  magnetic 
tape  units  Ml  and  M?  respectively. 

Phase  2 

After  the  remaining  IT's  are  processed,  Ml  is  rewound  and 
again  only  half  of  the  AITS  are  now  read  and  stored.  In  this  phase 
the  moving  sample  sum  for  every  224  samples  is  computed  and  if  less 
than  112  then  the  sum  is  complemented.  If  complementation  takes  place, 
in  effect,  the  A is  the  largest  value  and  retained  to  represent  that  zone. 
The  moving  sample  sum  is  denoted  as  ANUMA  and  can  be  expressed  as; 

224 

ANUMA(l)  = v AIT. 

i=l  1 
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ANUMA(i)  = -AIT(i)  + ANUMA(i-l)  + AIT(i+223) 

Every  224  samples  the  ANTJMA's  are  written  on  following 
the  BAITS.  The  remaining  AIT's  are  processed  and  recorded.  The 
BAITS  are  processed  identically  except  the  ANUMS  are  recorded  on  Ml. 

Phase  3 

In  the  third  phase  the  maximum  value  between  ANUMA  and 
ANUMB  is  chosen  to  represent  that  sample.  The  data  is  read  from  Ml 
and  M2  and  the  output  denoted  as  ANUM  is  written  on  M3. 

Phase  4 

The  fourth  phase  performs  three  functions: 

1.  Performs  peak  detection 

2.  Computes  the  statistics 

3.  Prints  results 

Range  peak  detection  is  accomplished  by  reading  in  from 
M3  through  ANUM  and  choosing  the  largest  value  every  672  samples. 

At  present  each  file  of  data  contains  33  signal  periods  $ therefore 
there  will  be  33  values  chosen.  The  computation  of  statistics  is 
accomplished  by  obtaining  a distribution  of  the  ANUM's;  1 . e . , it  will 
count  the  number  of  ANUMS  that  fall  into  113  theoretical  bins  numbered 
112-224  and  denoted  as  ISTAT.  Each  bin  represents  a threshold  value. 
Also  it  computes  the  percentage  of  total  that  an  ANUM  falls  within  each 
bin  and  all  lower  bins,  i.e., 

32 : 

Percentage  ($)  = v,  ISTAT  (i) 
i-1 


The  above  computation,  in  essence,  is  cumulative  distribution  of 

the  ANUM's  , and  the  results  shown  in  Figure  5 through  l4  are  a pictorial 
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representation  of  the  printed  results.  Flow  charts  are  shown  in 
Appendix  B.  The  program  in  its  present  form  utilizes  approximately 
54000  OCTAL  cells  including  library  subroutines. 

Inputs/Outputs 

1.  IT's  (mag  tape)  22400  or  22624 

2.  LSCALE  (punched  card)  specifies  input  S/N  ratio 

3.  IPEAK  (punched  card)  specifies  range  peak  detection 

a.  1 no  peak  detection  - print  distribution  only 

b.  2 no  peak  detection  - print  distribution 

c.  3 peak  detection  in  range  (every  672  samples)  and 
print  distribution 
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VII.  CONCLUSIONS 
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1.  The  primary  objective  to  adapt  the  WPP  simulation  program 
to  NUWC  facilities  was  accomplished. 

2.  In  order  to  obtain  a clearer  comparison  between  NUWC  and 
Sperry's  results  it  was  decided  to  triple  the  input  data  to  obtain 
greater  resolution  and  smoother  curves.  The  results  obtained  with 

the  adapted  program  are  nearly  identical  to  the  results  Sperry  obtained. 
Figures  3 through  l4  show  the  results  for  input  signal  to  noise  ratios 
(DB)  of  -®,  -6,  -4,  -3..  -2,  -1,  0.,  +1,  +2,  +3/  +4,.  +6,  respectively, 
for  both  the  adapted  and  the  original  program. 

3.  Presently  the  program  is  being  used  to  investigate 
processing  gain,  and  modification  to  the  program  is  underway  to 
investigate  range  binning  effects  on  processing  gain.  Plans  are  being 
made  to  modify  the  program  to  investigate  (l)  Mutual  Ship  Interference, 
(2)  Structured  Echo  Analysis. 
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FIGURE  5:  PROBABILITY  OF  EXCEEDING  THRESHOLD  VERSUS  THRESHOLD 


NUWC 

Sperry  Rand 


CONFIDENTIAL 


PERCENT  OF  F.ULL  SCALE 


THRESHOLD 


FIGURE  6;  PROBABILITY  OF  EXCEEDING  THRESHOLD  VERSUS  THRESHOLD 
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FIGURE  12:  PROBABILITY  OF  EXCEEDING  THRESHOLD  VERSUS  THRESHOLD 
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3.  Store  the  result  as  the  next  starting  point. 

4.  Convert  the  result  to  floating  point. 

The  initial  starting  value  is  a fixed  point  1. 


Notes:  1.  The  du may  argument  in  RANDM1  and  IRAND3  is  necessary  but 
its  value  and  type  are  unimportant. 

2.  A CS-1  output  324  with  the  relocatable  object  program  on 
cards  is  inserted  in  the  Fortran  deck  just  before  the 
$DATA  card. 

3.  The  routine  occupies  55q  locations. 
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i 


W'Vlil  i i i 

RANDM1 

RANDM2 

ENTRANCE 


INCREMENT 

EXIT 

ADDRESS 


COMPUTE  RANDOM 
INTEGER 

NEW  GUESS  = OLD 
GUESS  x 312  5 t 

536,870,909,  n 


ENTRANCE 


ENTER 

FUNCTION 

ARGUMENT 


STORE  INTEGER 
ARGUMENT  AS 
NEW  STARTING 
)INT  OF  GENERATO. 


CONVERT  INTEGER 
TO  FLOATING 
POINT 


ENTRANCE 


ENTER  CURRENT 
INTEGER  STARTING 
POINT  OF 
GENERATOR 
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GAUSS  - Pseudorandom  Number  Generator,  Normal  Density 

This  function  subprogram  generates  random  numbers  whose  density 
function  is  gaussian  with  zero  mean  and  unit  variance.  The  method  used 
is  a Mueller  Box: 

Given  two  independent  random  numbers,  and  , uniformly 
distributed  on  the  interval  0-1,  two  independent  gaussian-distributed 
random  numbers  are  obtained  via  the  following  algorithm: 
y1  = J-2  In  x^  cos  2 n x2 


y2  = J- 2 In  x1  sin  2 n x£ 


This  algorithm  was  obtained  from  Memorandum  9,  Random  Noise 
Generation  Routines,  February  19<o7>  prepared  by  E.  C.  Fraser, 

Stanford  Research  Institute,  Project  ESU  5830. 

The  method  is  the  same  as  that  used  in  the  1107  library  routine. 


given  by: 


fx  w = jrz  exp  (-x2/2> 


The  new  GAUSS  routine  is  written  in  Fortran  IV  with  the  calling 

sequence: 

R = GAUSS  (Dummy) 


Notes;  1.  The  routine  uses  the  uniform  random  number  generator  RANDM1 
twice  for  each  two  gaussian  numbers  requested. 

2.  The  routine  calls  the  SQRT,  ABS.,  ALOG,  SIN.  and  COS  routines. 
The  routine  occupies  102Q  locations. 
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COMPUTE 

[INTERMEDIATE  VALUE : 
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COMPUTE  SECOND 
ORDINATE  OF 
NORMAL  CURVE 
GAUSS  = 
C0S(2tt  xR  )xT 


/ COMPOTE \ 
SECOND  RANDOM 
NUMBER,  Rg 

'v  RAIIDM1  7 


COMPUTE  FIRST 
ORDINATE  OF 
NORMAL  CURVE: 
GAUSS  = 
SIN(?ttxR)xT 
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WAVE  PERIOD  PROCESSOR  NOISE  GENERATION 
FUNCTIONAL  FLOWCHART 
CONTROL  ROUTINE:  KAVEE 


NOISE  GENERATIOi 


INITIALIZE 
ACCUMULATIVE 
ARRAYS  TO  ZERO 


COMPUTE  & PRINT 
TOTAL  MEAN. 
MEAN  OF  SQUARES, 
VARIANCE.  AND 
DEVIATION  FOR 
ALL  BLOCKS  COMBINED 


COMKME  MOD- 5 
STORAGE  INDEX 


PRINT  CURRENT 
COMPUTATIONAL 
PARA^TERS 


'■'ARK  END  CF 
OUTPUT  ON  2 
TAPES 


/ GENERATE  1 \ 
BIiOCK  OF  noise' 
(16128  SAMPLES) 


EXIT  TO  > 
MONITOR  CONTROT, 
v SYSTEM  y 


COMTUTE  & PRINT 
VARIANCE  AND 
DEVIATION  FOR  1 
BIiOCK  OF  NOISE 


ACCUMUIATE 
MEAN  & MEAN  OF 
SQUARE  VALUES 
FOR  All,  BLOCKS 


COMPLETED 
101  BLOCKS 
OF  NOISE 


FUNCTIONAL  FLOWCHART 
COMPUTATIONAL  ROUTINE;  BLOCK 


NOISE  COMPUTATION 
START  


READ  & PRINT  11 
PARAMETERS  FROM 
CARDS:  SAMPLING 
TIME . BANDWIDTH 
& Z-TRANSFORM 
COEFFICIENTS 
OF  BUTTERWORTH 
FILTER 


INITIALIZE 
INDEX  VALUES 


COMPUTE  NEW 
SINE  AND  COSINE 
CHANNEL  VALUES 
(SEE  NOTES  FOR 
EQUATIONS) 


COMPUTE  5 NEW 
NOISE  SAMPLES 
FROM  SINE  AND 
COSINE  CHANNEL 
VALUES  AND 
LAGRANGE  CONSTANTS 


COMPUTE  LAGRANGE 
INTERPOLATION 
CONSTANTS 


MOD- 5 
INDEX  = 0 


< COMPUTED  \ 

AT  I EAST  \ 
16128  ) 
SA'tPLES  / 

V YES 

WRITE  1 BLOCK 
OF  NOISE  ON  EACH 
OF  2 OUTPUT  TAPES 


TRANSFER  EXTRA 
VALUES  FROM  END 
OF  OLD  BLOCK  TO 
BEGINNING  OF  NEW 
BLOCK  OF  NOE  E 
SO  THAT  THE  NOISE 
FUNCTION  IS  CONTINUOUS 


COMPUTE  MEAN 
OF  SUM  & MEAN  OF 
SUM  OF  SQUARES 
FOR  CURRENT 
BLOCK  OF  NOISE 
(16128  SAMPLES) 


ADJUST  TJffi  MAJOR 
PROCESSING  INDEX 
ACCORDING  TO  HOW 
MANY  SAMPLES  ARE 
TRANSFERRED 


exit  to 

CONTROL  ROUTINE 


Kj  Jl'«<  ’ 


55 


. 1 »«- 


WAVE  PERIOD  PROCESSOR  NOISE  GENERATION 
PROGRAM  FLOWCHART 
CONTROL  ROUTINE:  KAVEE 


)ISE  GENERATION 
START > 


INITLALIZE 

DATA 


CLEAR  2 5 WORD 
DOUBLE  ARRAYS 


XS,  XC 


CLEAR  2 3 WORD 
REAL  ARRAYS 


YS,  YC 


CLEAR  2 REAL 
VARIABLES 

SUMZ,  SSQZ 


EJECT  PAGE 
ON  PRINTER 


J = REMAINDER 


2x  ( NBIiOCK- 
5 


MOD. 5 STORAGE  INDEX 

y 

PRINT  CURRENT 
COMPUTATIONAL 
PARAMETERS 


' PERFORM 
NOISE 

COMPUTATIONS 
, BLOCK 


COMPUTE  VARIANCE 
& STANDARD 
DEVIATION  FOR  ONE 
NOISE  BLOCK 
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V  
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NOISE  BLOCK 
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NOISE  BLOCKS 
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V 


COMPUTE  TOTAL 
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1 
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PROGRAM  LISTING: 

Noise  Subroutine 
Single  Precision 


p r'  ■;  ’Tm 

uu i ti  ijl i j I InL 

PAIR  SYSTcvl  ANALYSIS  wave  PERI0D  Pr0CEsS9r 
FILTERED  3AJSSIAN  N0ISE  - FIVE  P0LE  8UTTERW8RTH  FILTER 
O0UBLE  PRECISION  XS(5).XC(5) 

0ImENSI3N  YS(3).YC(3) 

INTEGER  kMDSET 

data  NTRT/o/.RNDSET/01/.LNgTH/288/»MULTP/56/.NBLKS/1O1/»BLKS/1O1./ 


600  F0RMAT ( 1 nl //) 

601  F9RMAK5UX.  19).  8X  , 6 1 2/2  ( 5D23.  1 6/)  . 2 ( 3E2  3 . 8/)  ) 


602  F0RMAT(4o24. 3/) 

603  F0RMAT<5A, I 5.8X.912,8X,012) 
PAUSE 

REWIND  1 


REWIND  2 
03  5 Is*, 5 

X5{  I )=0.oDO 
5 XC  ( I ) =0.  oDO 


03  10  1=*, 3 
YS ( I ) =0, u 
10  YC(I)s0.u 
SUmZ=0.0 
SSqZs0*  0 
WRlTE(6,u00) 

NBL®Ck  = 10* 

09  15  NBl3CX=1»NSLKS 
J s M8DU*(NSL9CK-l>.5)+0 

WRITE (6,oqi )N0L9CX.LNGTH,MULTP. J, N TR Y. RNDSET. XS, XC. YS, YC 
PAUSE 

CALL  BL9C<(N9L3ck.LnGTh. J.NTRY.RNDSET.XS. XC* YS, YC, YBAR, YMSQ) 

VARY=YV5j-Y3AR*#2 

Sl3Y=SQR« (VARY) 

WRITE (6, 032) YBAR, YMSQ. VAR Y, S I GY 
SUmZ=SUM4+Y3AR 


ss3Z=sso4  + rvisa 

15  CSsjTInUE 

ZBaR  = S'J\U/3L<5 

ZM50  = SSO*./?L<s 
VARZ=ZVSo-Z3AR*#2 
SIgZ=S2RI (VARZ) 

WRlTE(6,032)Z3AR,ZMsQ,VARZ,SlGZ 

WRITE ( 5,031 )N3L3CK,LNGTH, MULTP.J.N TRY. RNDSET,XS,XC,YS,YC 

end  file  * 
rewind  1 
ENc  File  2 
R E a 1 \ D 2 


S T ; P 
END 
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CGm’iULu  I ir\L 


Wmi 


SUd^aUTlNu  3L3CK(N8L0CK.LNGTH. Jf NTRT, RNDSET.XS.XCt  YS. YC. YBAR* YMSO) 
0 1 MENS  I a*  Li T3&) *L2(Sb) 

C FILTERED  GaJSSIaN  N0JSE  - FIVE  P9LE  8UT1ERW0RTH  FILTER 

D9UBLE  P*ECIS13N  XS ( 5 ) . XC ( 5 ) . TEMP. I , 3rt , C i . C2. C3. C4. C 5 
Dl MENS  IS*  YS(J>.  YC  (3)»Y(  16132)  *G  (209) 
integer  ksdset 

601  F0RMAT(5!2XI 18)) 

503  F0RMATC2UXE23.  16>/3XE23.  l6t2(4XE2P.  l6)/2(4XE23,  1 6) /4  (4XE  1 5*  8)  ) 

603  F0RMAT  < 4rt  T sE23.16.4H  BW=E23.16/ 

14h  C1=E2R. 15.4H  C2=£23.16.4H  c3=E23.l6/4H  C4sE23.l6.4H  C5se23.16/ 
24H  C6=E14,9.4H  C7=El4.8»4H  C8=E14.»»4H  C9=E14,8/> 

!F(NTRY)a,5.  IS 

5 REaD  (5.^33) T.BW, Cl, C2.C3.C4.C5.C6.C7.C3.C9 
WRlTE(6.o33)T»3.V,Cl.C2.C3.C4.C5.C6.C7.C3.C9 
MJlTP=55 
MLTM 1 =MJwf3"l 


c 


15 


20 

25 


NJmDR=MJuT?*LNGTH 
D£m0M  = N J143R 

TH9EE  PSiMT  LAGRaNGe  I NTERPSL aT 1 9N  C9NST ANTS 

SAls-.OS«. 95105652 

SA2=*. 96*. 95105652 

SA3=*. 12*. 95105652 

CAis-. 03*. 30901699 

CA2=+.95*. 30901699 

CA3=*. 12*. 30901699 

SBis-. 12*. 53^8525 

SB2=+.84*. 55773525 

SB3=+. 23*. 53778525 

C3 l =-. 12*. 53901599 

C32=+. 84*. 30901699 

C33s+. 25*. =0901599 

NTt?Y  =1 

IFlAG=0 

I^iAX  =3245 

IF(J.EQ.^J33  T 3 25 

DS  20  Isa.j 

Msi+NUVBK 

Y(!)=Y(Ma 

IF  ( J.GT.  4)  I v;ax=i*ax-i 
11=0 

D 9 3C  l=u, I VAX 

TEmP=c1*^5(5)-C2*XS(4)+C3*XS(3)-C4*XS(2)+C5 

XS(l)=X5i2) 

XS(2)=XSv3) 


*XS(1)+GAUSS(2) 


XS(3)=XSv<i) 

XS(4)=X3v5) 

X 5 ( 5 ) =TEm»  : i 

TEmP=c:*a!:(5)-C2*XC(4)*C3*XC<3>“C4*XC(2)+C5  «XC  ( 1 ) +GAUSS  ( 1 > 

XC(l)sXCl2)  v 

XC  (2)  =x:  v.3) 

XC  (3)  sXO  U) 

xc ( a ) =x: v 5 > 

XC(5)=T£rtD 
YS;  i)=Y5<.2> 

Y 5 { 2 ) = Y 5 v 3 ) 

YH  ;3'  =C?*-;'T3LiXS(4)  )+c7*SNGL(XS(3)  1+c8*5N3L<XS(2)  )+C9*SNGL(XS(l)  ) 


M=5* » ♦ J ♦ 2 

Y (n)  s-(5«»*rsc 1) +SA2*YSC2) + 5M*YS< 3) ) 

2 ♦CA3»YC(l>*CA2*YC(2)*CAi#YC(3) 

M=5* f ♦ J* J 

Y(«')=YCU) 

m=5*:+j+4 

Y ( »»)  s 5^i*rSU)+SA2#YS<2)+S^3*YS<J) 

4 ♦CAi#YC<i)*CA2#YC(2)+CAJ#YC(3) 
M=s# ' * J + 'j 

y(M)s  SDi*rSU)+S32«YS(2)+593*YS<3) 

5 -(CB;*YCC  .)-*“C32*YC(2)+CB3*YC  ( >)) 

30  C5^T:r;Ur 

WrtJTECn  (V(‘J) » M = l.NUMBt?) 

•I  ft  i T E C 2 ) (Y(M)»  M = i.NUMflR) 

* I w | 

roi  R -DSET=  : •< <\*i03  c ; > 

SwVY  = C*  2 
S5qY  = C*  2 

Otf  03  K =O.MUl?U 
D«  32  Kisi.L'-iGYH 

T it*  < 1 

32  Gl*i)=YU) 

D 5 a 3 I=i,-^sSIH 
Svv'Y  =5  JMY  *3  ( I > 

45  5S2Y“SS2Y^3 ( J ) **2 
85  C-VTiNUE 

Y2,rK-5,J,,Y/3EMaM 

Yv.SOsSS  jY/0£*Ja^ 

RETURN 
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Ui'I 

111  ilki  I 


r-v  — (•  i-j-j  0 1 

UiiuilAL 


r 


! 


I 

l 


100 


FUN,C7:,'‘J  liAjgsCAA) 
ir  (IFLAu  . N£.  C)  30TQ  100 
i Flag  = • 

tkanw  s ua'jd** (oummy) 

TcmP  = 5-^T(2.*ABS(aLQG(TRaNR))) 
ANGLE  s *AN3Mi  (DUMMY) 

GAUSS  = -iM(6.283185*ANGLt)*TEMP 

W£TU*N 

1 r L AG  s w 

GAUSS  = w3S(6»2a3lS5*ANGL£)«TEMP 

rtcTU^N 

end 


r-  'Cir'Ti'.'TjM 

V-* * _/ 1 « * i J I I/-.L 


ct 


I 


! i 
! 


I i 

• -j 


PR8GRAMS  Ls*050  AUD 
MAIn  i 

BLICK 

GAUSS 

* RANoM 
A I 6 2 

dsle 

CBS 
SIN 
E I 2 1 
E 1 22 
AL9S 

EXP 

ABS 

M30 

* SORT 

* SNGL 

* MClSl 

* MC  I 62 

* AS  I 22 

* M I 22 

* D I 22 

* N I 22 

* C I 2 1 

* C I 7i 

* AINT 

* I AT 

* A I 7 1 

* Cl  26 

* AS  I 66 

* M I 66 

* C I 62 

* C I 6 1 

* C I 12 

* ZI70 

* CIvT 

* ARLl 

* exit/ist 

* F3RTI/3 

* FIER 

* IPA 

* F3RT/MS-1 


?RE  LIVITs 


00600 

- 

01201 

01202 

- 

01276 

01277 

- 

03211 

03212 

m 

43561 

43562 

■> 

43651 

43652 

m 

43665 

43666 

- 

43726 

43727 

- 

43743 

43744 

- 

43745 

43746 

- 

43757 

43760 

- 

43774 

43775 

- 

43776 

43777 

- 

44120 

44121 

- 

44144 

44145 

- 

44 1 63 

441&4 

- 

44  i65 

44166 

- 

44276 

44277 

- 

44311 

44312 

- 

44407 

44410 

- 

44423 

44424 

• 

44534 

44535 

- 

44557 

44560 

- 

44573 

44574 

- 

44620 

44621 

- 

44621 

44622 

- 

44736 

44737 

m 

45005 

45006 

- 

45026 

45027 

• 

45055 

45056 

- 

45205 

45206 

- 

45277 

45300 

- 

45376 

45377 

- 

45403 

45404 

- 

45431 

45432 

- 

45440 

45441 

- 

45467 

45470 

- 

45525 

45526 

- 

45543 

45544 

- 

45605 

45606 

46034 

46035 

- 

46232 

46233 

- 

46266 

46267 

- 

46322 

46323 

. 

46367 

46370 

. 

46405 

46406 

- 

46461 

46462 

- 

46517 

46520 

- 

46521 

46522 

- 

52242 

52243 

- 

52327 

52330 

- 

52345 

52346 

- 

53003 

CONFIDENTIAL 
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COiNriuti'iiiAl 


— ■mWlM 


E'lTRY  POINT  AND  C«MMSN  BL0CK  LeCAT 1®^S 


NAME 

L0CATI9N 

fblank 

ooooo 

maini 

00602 

1 PA 

52330 

|RE 

52162 

| MR 

46612 

|ND 

47340 

mbd 

44576 

1 ix 

47276 

1 10 

47265 

1 IF 

47272 

BLBCK 

01301 

E 1 2 1 

44147 

S 1 22 

45056 

SORT 

44622 

A 122 

45064 

D 1 22 

45300 

1 EN 

52155 

1ST 

46520 

1 AT 

45470 

| RD 

46605 

Cl  12 

46323 

M 1 22 

45206 

gauss 

43564 

M 1 66 

46035 

S 1 66 

45606 

Al66 

45627 

A 1 62 

43731 

sngl 

44737 

1 2 1 

51251 

INE 

51330 

I RAnD3 

43722 

RANqMI 

43666 

ALBg 

44314 

aBS 

44562 

SlN 

44001 

ces 

43762 

RANdM2 

43712 

dble 

43750 

cl  26 

45544 

C 1 2 1 

45404 

N 1 22 

45377 

E 1 22 

44170 

EXP 

44426 

aint 

45441 

FlER 

52243 

C 1 7 1 

45432 

2 1 To 

46370 

A 1 7i 

45526 

gRRgR 

52255 

idint 

44754 

Cl62 

46233 

MCI  62 

45027 

C 1 6 1 

46267 

MC  1 6 1 

45006 

aRli 

46462 

ClVTARGLl 

46444 

AR6L1 

46505 

cl  VTAQ 

46424 

ClVTO 

46406 

I ifa 

45502 

■ -l..  :A«**  ...  — J Ml 

oGiir'IDENTIAL 
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M«N|  | 

74506 

iba 

52150 

I IE 

47260 

|0E 

4652F 

1 EC 

46522 

|EF 

46566 

ibreti 

46765 

ibtbli 

52636 

WRITE* 

52364 

SETeSF 

52346 

readi 

52361 

BACKSPACE! 

52353 

endfilei 

52356 

rewinoi 

52367 

E 

PROGRAM  LISTING: 

Noise  Subroutine 
Double  Precision 


► 

\ 


CONFIDENTIAL 

PAIR  SYSTtM  ANALYSIS  WAVE  PERJ8D  PR8CESS9R 
c FILTERED  3AJSSUN  N0ISE  - FIVE  P8Lt  8UTTERW6RTH  FILTER 
IMPLICIT  D»JBLE  PREcISISn  (B» S» J-Z) 

08UBLE  PK£  C I S l 9N  XS(5).XC(5) 

D I MENS  1 3>»  YS<J).YC(3> 
integer  *yoset 

DATA  NTRl/3/»RNDSET/ei/»LNGTH/2B8/»MULTP/56/»NBLKS/10i/ 

DATA  BL<a/lDl»ODOO/ 

600  F 9RM  A T < 1 rtl //) 

601  F0rmAt(5A1X.  19).  0X.6l2/2(lX.5D23.1t>/).2(oE23,8/)) 

602  F3pvAT<4»*24.B/> 

603  F9rmAT<5a, I 5» 8X, 912. 8X, 812) 

PAtjSE 

REMIND  1 
REMIND  2 
09  5 1*4,5 

X5(  1 >=0.000 
5 XC(I)sO.w0.0 
03  10  1*4,3 
YS(J)sO. j 
10  YC(I)s0.o 
SUyZ*0.0 
SSqZ=0.0 
ARlTEtb.oOO) 

503  F8rmat(3a, ?12/3<3X,d23. 16)/2(3X.D23. 16)/3(3X.D23.  16)/2(3X.D23, tr) 
1/3<3X,02.>.  *)/3< 3X.  023.0 >/> 

REaD(5*Sw3)RNDSET.(xS(|). 1*1.5).  ( XC ( l ) • I = 1 . 5)  • (YS(l)» 1*1,3), 

1 c Vc  < I ) • 1 = 1,3) 

c NBlSCk  = 10* 

D9  ]5  N3*.3C<s)»NBLKS 
J r V0D(C.('4BL8C<-1).5)A'O 

WRjTE(6.  tO! )NBLSCX,lNGTH,MULTP. J.NIRY.RNOSET.XS. XC, YS, YC 
CALL  BL3^< (NBLSCX.LNGTH, J.NTRY.RNDStT. XS.XC. YS. YC.YBAR.YMSo) 
VArY=YYS--Y3AR»*2 
SI5Y  = D*3R  T (VARY) 

WRITE (6.^32) Y3AR, YMSQ, VARY, S I0T 
SU\fZ*SUvU*Y3AR 
SSoZ  = S^04.*Y  MS3 
15  C3\TINU£ 

ZBaR  = SUV»4./3l*5 

ZmsQ=SS3a/?lXS 
VArZ=ZMSm-Z3AK..2 
SIgZ  * o-ormtarz) 

WRITE .?.w32)Z6AR,ZMS0. VARZ, SI GZ 

WRITE (6.u01)NBL8CK,LNGTH.MULTP. J.NIRY.RNDSET.XS. XC, YS, YC 

EN0  FILE  I 

REhImd  1 

END  FILE  2 

REWIND2 

step 

ENq 
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601 

503 

603 


15 


20 


25 


1NTEHP9LAT10N  CONSTANTS 


S08R9UT 1 N&  8L»CK(NBL0CK.LNGTH.  J»  NTR Y,  RNOSeT. XS, XC. YS,  YC. YBARt  YmSO) 
IMPLICIT  0?J3LE  PREcISIBN  (A-H.0-X.4) 

OImENSIB.'*  lH56).L2(55) 

FIlTERED  3 * JSS I AN  N0ISE  - FIVE  P8LE  9UTTERN0RTH  FILTER 
OeuflLE  PKECISI3N  XS ( 5) . XC < 5) # TEMP, 1 . 8*. C 1 . C2. C3. C4. C5 
DOUBLE  PrtECiSIS.N  YS(3>,  YC(3)»  Y BAR*  YMSQ 
01m£NSI9^  Y(16132).  G(288) 
integer  rndset 

F8pvlAT<5l2Xl  18)  ) 

F3rv<At<2UxE2J.  161/3XE23.  l6,2(4XE2J.  l6)/2(4XE23.  16) /4  (4XE 1 5*  8)  ) 

F A T ( 4rt  T =E23, 16, 4H  BWsE23.16/ 

14h  C1=E2«>,  16. 4H  c2  = e23.16,4H  C3*E2>>.  16/4H  C4*E23.16.4h  C5se23.16/ 
24H  CfisE 1 p, A*  4H  C7SE i 5. 8*  4H  C8=El5.8*4H  C9=El5.8/> 

IF(NTRY) j,5. 18 

REaD  (5.333) T.BH, Cl, C2.C3.C4.C5.C6.C7.C8.C9 
WRITE (6,033) T . Bw, C 1 , C2. C 3. C4. C 5. C5 . C 7. C d . C9 
MUlTP=56 
MLTM 1 =MUuTs-l 
NUmBR  = MJ»-To»LNGTH 
0En8M=NJ43R 

Three  P3*  nt  lagrange 

SAi=-, OB*. 95105652 
SA2=*.96». 95105652 
SA3=+. 12*. 95105652 
CA l 08*. 30901699 
CA2a*,96#, 30901699 
CA3r+, 12*. 30901699 
SBia-,i2».5377B525 
SB2=*. 84*. 53778525 
SB3=*. 28*. 53778525 
Cd l®-, 12*. 83901699 
CB2=*. 84*. 30901699 
CB3S+. 28*. 33901699 
CALL  RANw'(2  (RNDSET ) 

NT9Y  si 
IFlAGsO 
ImaX  s32*5 
1F(J.eQ.M)33  T 8 25 
03  20  l=i, J 
Ms  | ♦MjmBk 
Y<l)=y(vf) 

IF(J.gT.O  IMAXsimAX-1 
I 1=0 

03  30  I s j, IMAX 

TEmP=C1*as(5)-C2*XS(4)*C3*XS(3)-C4»XS(2)+C5 
XS(1)=XS12) 

X5(2)sXSl3) 

XS(3)sXSU> 

XS(4)rX3(5) 

X5  ( 5 > = 1 

TEmP=C1*a:(5)-C2*XC(4)*C3*XC(3)-C4*XC(2)+C5  *XC ( 1 ) ♦GAUSS < 1 ) 

XC(l )sXC(2> 

XC(2>=x:c3) 

xc(3)sx:w) 

xc (4) sXC  («> 

<C  (5)  =7Ef13 
Y 5 ( l ) = Y 5 V 2 } 

Y5(2)=Y5C5i 

)S(35  =C‘-*\S(4>*C7*Xs<3)«-CQ#XS(2)  ♦C9*XS()) 


•XS(1)*GAUSS(2) 


yc(Dby:;2) 


I 


I 


. ] 


yc(2;=y:vj'. 

YC  (3  5 =C6.»\C<4>*C7*XC(3>*C8«XC(2)4-C9*XC(  ;) 


-IQQl 


CONFIDENTIAL 
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h»5*  J 

Y<M)a-CS«>3*YSU)4S92*YS(2)*S8l»YS(J) 

| ♦C83»YC(i}>CB2#rC(2)*Cdl*YC(3)) 

M»5« • 

Y (.<)=- <Sa3.YSU>*SA2»YS< 2 >+SM*YS<3>  ) 

2 ♦CA3»YC(l)*CA2*YC(2)*CAl.YC(3) 
m*5* \ +J+4 

Y<*0=YCU) 

M*5*  [ ♦ J*4 

Y (m)  a S«1*YSU)  ♦SA2#YS(2)+SA3«YS(  J) 

4 ♦CAi»YC(i)*CA2*yC(2)*CA3#YC(3) 

M»5* I+J+S 

Y(m)b  5t»l*YSU)*S32*YS<2)+S33*Y5(3) 

3 -tCBi*yc<;>*CB2*YC(;.»>*C93*YC<J>> 
30  c^vt;noe 

- i.MUMBR) 

= It  NUMBr) 


cor;?iCEWTi,u 


(YCK). 


TQl 


32 


I Tf:  C i 3 
rtSjTE(H) 

X 1 s I I ♦ X 

R.«iDStT=I'<A*J03U> 
SjvY"C«  3 
SS3Y=CO 

DS  95  =0.v|LTM1 

03  32  KisuLNGTH 
H*K*L.\GTri*ia 

G(ft4>sY  (,,) 

OS  45  isi.L'JGTH 
S w Y = 3 ’J  '1 Y * 3 ( I > 

45  SSoY3Si>3Y  + GU3**2 


CONFIDENTIAL 


FUNcTISN  «A'j5StAA> 

implicit  dsjble  precisisn  <a»g» t> 
O9(J0LE  PK£CISI9N  GAUSS 
IF  ( IfLA*  .SE*  0)  G8T0  100 
IFLAG  « A 

TRaNR  * kandmi (DUMMY) 

TEvP  s Do33T<2.*DA9s<DL9G(TRAnR>>> 
angle  s WaNONI  (DUMMY) 

GAUSS  * w3IN(6.293185*ANGLE)*TEMP 
RETURN 

I F lA  3 * w 

GA'JSS  = wC?5<6.283185*ANGLE)  *TEMP 

return 

EN0 


PPOGRAmS  igAocf}  A N Q 

UNITS 

MAInI 

00600 

m 

01347 

01350 

m 

04471 

BL8CK 

01472 

m 

03576 

03577 

m 

44627 

GAUSS 

44630 

m 

44735 

44736 

m 

44757 

• randm 

44760 

m 

45020 

E 1 6 1 

46021 

m 

45042 

45043 

45046 

El  66 

45047 

m 

45214 

45215 

m 

45234 

DC9S 

45235 

m 

45253 

45254 

- 

45257 

OSIN 

45260 

- 

45424 

45425 

• 

45475 

DL9G 

45476 

- 

45612 

45613 

45661 

DEXp 

45662 

- 

46031 

46032 

• 

46075 

DM9o 

46076 

- 

46125 

46126 

- 

46127 

DABS 

46130 

- 

46144 

M 1 62 

46145 

• 

46161 

46162 

- 

46163 

03LE 

46164 

• 

46175 

MiO 

46176 

- 

46222 

46223 

• 

46223 

• DSOrT 

46224 

- 

46327 

• M 1 22 

46330 

• 

46421 

• C » S 1 

46422 

- 

46431 

• DINT 

46432 

- 

46467 

• 1 AT 

46470 

• 

46525 

• C 1 26 

46526 

- 

46567 

* ZI8g 

46570 

• 

46614 

• A | 8 1 

46615 

- 

46645 

* AS  I 66 

46646 

- 

47074 

• M 1 66 

47075 

- 

47272 

• 0166 

47273 

- 

47500 

* N 1 66 

47501 

- 

47511 

* C 1 62 

47512 

- 

47545 

* C 1 6 1 

47546 

- 

47601 

* C I ft  6 

47602 

• 

47627 

* CIVT 

47630 

- 

47703 

• ARli 

47704 

. 

47741 

• EXJT/IST 

47742 

- 

47743 

• F3RTI/8 

47744 

- 

53464 

* FIER 

53465 

• 

53531 

• IPA 

53552 

- 

53567 

• F8RT/MS-1 

53570 

- 

54225 

COfJFiDEw  l in 


i 

nL 


EMRY  P3INT  AND  CBMmS'J  BL9CK  LOCATIONS 


NAME 

LfiC AT  1 9N 

FBLaNK 

ooooo 

MAlNl 

00602 

|PA 

53552 

|RE 

53404 

Cl  26 

46529 

| NR 

50034 

| NO 

50562 

| RO 

50027 

1 lx 

50520 

HO 

50507 

490 

46200 

818CK 

01474 

E 1 9 1 

45023 

S 1 ®6 

46646 

DSQrT 

46224 

A 1 9ft 

46667 

0196 

47273 

|EN 

53377 

1ST 

47742 

1 AT 

46470 

Cl  16 

47602 

4 122 

46330 

RANDM2 

45004 

gauss 

44632 

M 1 96 

47075 

C 192 

47512 

1 21 

52473 

|NE 

52552 

jRA<gD3 

45014 

RANdMI 

44760 

0L8g 

45500 

OABS 

46132 

M 1 ®2 

46147 

OSIn 

45262 

oces 

45237 

E 1 96 

45051 

0E*P 

45664 

0M90 

46100 

FIER 

53465 

C 1 9 1 

47546 

N 1 9& 

47501 

C 1 6 | 

46422 

ZI80 

46570 

a 1 a 1 

466(5 

dint 

46432 

OBLe 

46166 

eRRpR 

53477 

ARLi 

47704 

ClVTARGLl  47666 
AR6Ll  47727 
ClVTAQ  47646 
I I F a 46502 
ClVTO  47630 
m9N||  74506 


-■ 


*■  - ■ 

i*«et 1 

5020? 

CONFIDENTIAL  1 

I ®tbl i 

54060 

«RlTE 1 

53606 

SETeSF 

53570 

ftf  AOI 

53603 

1 

BACKSPACE  » 

53575 

■ 

ENDFUEI 

53600 

• 

REk | NO | 

53611 

.... 

CONFIDENTIAL 


f 


PROGRAM  LISTING: 


Wave  Subroutine 
Single  Precision 


i 


1 


C NAVE  PeRISO  MEaSORMENTS  sf 
DI  MENS ! ® TT 
DImENSIS*  IT(224#12) 

DImENSIS.'.  SPNU6128) 

EQUIVALENCE  Uf<l).TT(l>> 
data  L»UI/224/.  LSCaLE/12/.  FnP/12.0/ 

Data  P 1 /9.  I41t>927/.  F 0/20000. 0/. FBW/4 40. Q/. DELT/O. OOOOIO/. NP/ 12/ 
DATA  NBL<N5/10l/,LNGTH/283/.MULTP/5t>/ 

DATA  RMS/ 1. 0004576/. SORT  2/1. 4 1 42 13&/.TP/0. 16128/ 

DATA  8NEU3/1. 12201 S/.THRD8/ 1.4 1253//. TENDd/3. 1622777/ 

REWIND  l 
REWIND  2 
REwlNO  3 
PAUSE 

DBlPHF>U  122018*«l.s 
Tn§DBs1.A2201B**2 
F3dB*1.  U2318»*3 
F4DB*TW8J3»»2 
F5dB*1.  U20iB»«*5 
F6DB  * T«^D3»*2 
F 1 2DBsF  6“3  * »2 
Fl8DR=F6w3«*3 
F2aDBsF 1*D8**2 
F3oDBsTE‘'*33»#3 
611  F8RMAT(5i2Xl 18)/) 

c change  thu  for  change  in  lbut 

c S/.M  5 peak  »I3NAL/SQRT(2.0)/RmS  NBISE 

C S/N  RAT  I 9a 

SCaLE<1)=0. 0 
SCaLE(2>s^\»5*S0RT2/F6DB 
SCAL£(3)a^v(3»SQRT2/F4DB 
SCAlE(4)*:?vt5*S3RT2/F3DB 
SCALE  (S)*Rv(S*S0PT2/TW3DB 
SCALE(6)i^MS*!»0RT2/eNtDB 
SCaLE(7)»^mS*SQRT2 
SCALE  <8)s}M3«*SQRT2*eNEDB 
SCALE  (9)s^vl3«SQpT2*TW0DB 
SCALE(10>=RMS«SQRT2*F30B 
SCalE(11)=Rv(S*SQRT2»F4D8 
SCALE  ( 12>=Rv(S*SQRT2*F6DB 
MLTM1  sviJLlP-l 

lit  *LSCALE»L®uT 

TW0  P|  S2.0*P1 

Fl  *3.5*FBW/Tp 


signal  plus  nsise 

(224.  12).  SC ALE  1 1 2 ) 


CONFIDENTIAL 


NUMBR  *mJL1P»LNGTH 

c Change  thu  fjr  change  in  l9ut 

NUm3R  * A 5 1 28 
FNuMBR  sNJMBR 

TNUMBR  s(FNUMBR41.0)/2.0 
• SQL  * 34.0*FO»FO/(FNP*FB.1) 

c DIGITAL  wLSCK  FREQUENCY  2.424  MEGACYCLES 
TSCALE  sDELT*8SCL 

L2  *LTT 

L6  sLTT 

C8  60  N*4»  V 3 L K S 
READd)  »3N 
REWIND  2 
NR|TE(2)  SPN 
D»  55  L**,LSCALE 
lF(L»EQ.4)3t  T9  23 
REWIND  2 
REaD<2>  ®PN 
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CONFIDENTIAL 


i 

•j 


CONFIDENTIAL 

c 19.2  KVD*  RANGE  (CHaNGE  FSLLSwING  CARD  F9R  DIFF  RB  BIN) 
IF(Meo(N*3)*NE.l)  G9  T9  30 
23  09  25  Jst.NJUdR 
SJ  = J 

TX  s DEL»*(SJ-TNUMBr) 

c generaat * sn  9F  linear  fm  pulse 

s s C®S(  » A'9PI*TX«(F0*F1*TX)  ) 

25  SPn(J)  a 53N(J)+SCAlE(L)«S 
30  CONTINUE 

Da  55  |s*,l30T 

C CiMpUTE  STAR!  = FIRST  P0SITIVE  ZER0  CROSSING  rtITHiN  SAMPLE  INTERVAL 
J=72*J-7* 

40  1F(SPN(J)  >41.42*43 

41  J = J*  l 

IF(SPN<J>  J4l.45.46 

42  J*J*l 

IF(SPN(J>  )41. 45*44 

43  J=J> 1 

GO  T 3 40 

44  Jaj-1 

45  STaRTsJ 
68  T®  50 

46  P LlNrSPN(J>/(SPN(J>-SPN(J-l ) ) 

SJ  s J 

STaRT  = (UP  LIN-1. 5>*P  LIN  + 0.5)«U. 566-1. 0)*P  LlN*SJ 
C CSJnT  NP  POSITIVE  ZER®  cR®S5JNGS  aNO  THEN  C8MPUTE  ELAPSED  TIME  TT(I.L) 

50  D8  53  1CH3SS=1*NP 

51  Jaj+1 

IF(SPN(J)  >52. 52*5i 

52  J = J ♦ l 

IF(SPNCJ>  >52.53*53 

53  CSnTInUE 
S J s j 

TT(I.L>  a SJ-START 
IF(SPN(J1  >54,55*54 

54  P LlNsSPt'i(J>/’(SPN(J)-SPN(J-l)> 

TT(t.L)  =<<<P  LlN-i.5)*P  LlN+0.5>»0. 566-1. 0>*P  LIN+TT(1,L> 

55  CONTINUE 

C MSoULS  3c  C3UNTER  0UTPUT  WITH  SLIDING  PRESET 
09  64  JKS 1 , LSCALE 
I SET  -0 

D©  64  I s l » L®UT 

c Change  thij  fsr  change  in  l®ut 

IF  <M0D(i,mO  ,E0.  0>  ISET  a ISET+1 

54  IT  (l.JK)  s M0D(INT(TSCALE*TT(I, J<> >*lSET,32>+0 
NR  I TE ( 3 > IT 
NRlTE<6,oii>  N 

pause 

60  CONTINUE 
REWIND  l 
END  FILE  3 
REWIND  3 
ST©P 
ENq 
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P98GRAMS  L0ADCD  AND 

C9RE  LIMITS 

MAI  N 1 

00600 

- 

01604 

01605 

• 

46556 

css 

46557 

m 

46573 

46574 

m 

46575 

SIN 

46576 

- 

46717 

46720 

m 

46743 

E I2l 

46744 

m 

46762 

46763 

- 

46764 

El  22 

46765 

- 

47075 

47076 

- 

47110 

alsg 

47111 

- 

47206 

47207 

- 

47222 

EXP 

47223 

- 

47333 

47334 

- 

47356 

ABS 

47357 

- 

47372 

M3D 

47373 

47417 

47420 

- 

47420 

1 NT 

47421 

47432 

* AS  1 22 

47433 

- 

47562 

• M 1 22 

47563 

- 

47654 

• D 1 22 

47655 

m 

47753 

* N 1 22 

47754 

- 

47760 

* CI21 

47761 

m 

50006 

• C 1 71 

50007 

m 

50015 

• A I NT 

50016 

- 

50044 

• 1 AT 

50045 

- 

50102 

* A 1 7 1 

50103 

50120 

* C 1 12 

50121 

- 

50165 

* Z 1 70 

50166 

m 

50203 

• CIVT 

50204 

- 

50257 

• ARLl 

50260 

m 

50315 

• exit/ist 

50316 

m 

50317 

• F3RTI/S 

50320 

m 

54040 

• FIER 

54041 

m 

54125 

* 1 PA 

54126 

*• 

54143 

* F3RT/MS-1 

54144 

- 

54601 

1«ft 


ENTRY  P9INT  AND  C®MM9N  BLSqK  LBCaTISNS 


SAME 

LBCATIBN 

PBLaNK 

00000 

M A I N 1 

00602 

I RE 

53760 

|PA 

54126 

E 1 22 

46767 

£l2l 

46746 

M 1 22 

47563 

0122 

47655 

C>  12 

50121 

A 122 

47441 

|RD 

50403 

1 21 

53047 

|NE 

53126 

| NR 

50410 

mSD 

47375 

S 1 22 

47433 

C6S 

46561 

INT 

47423 

1 1 x 

51074 

|ND 

51136 

1 en 

53753 

|ST 

50316 

1 AT 

50045 

SIN 

46600 

ABS 

47361 

C 1 2 t 

47761 

N 1 22 

47754 

ALBg 

47113 

EXP 

47225 

aint 

50016 

PIER 

54041 

cl  Tl 

50007 

ZlTO 

50166 

AlTl 

50103 

ARLj 

50260 

Cl VTARGL1 

50242 

ARGlI 

50303 

cl VtAO 

50222 

ERRsR 

54053 

civtq 

50204 

i'Fa 

50057 

M0N  | 1 

74506 

|BA 

53746 

1 IF 

51070 

1 ID 

51063 

1 IE 

51056 

IDE 

50325 

|EC 

50320 

IEF 

50364 

10RETI 

50563 

19T0LI 

54434 

WRITE  1 

54162 

SETE0F 

54144 

readi 

54157 

BACKSPACE! 

54151 

endfilei 

54154 

re*indi 

54165 

Uui  i i lULi  <1  rlL 


i 
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por, 


Lo  i 


nave  period  *eaSurments  $f  signal  plus  nsise 
implicit  dsjble  precision  <a-h.0-r»u-z> 

DOUBLE  PRECISION  TP. THRDB. TENDS. TWtfOB* TW3P I • TNUMBR. TSCALE . T* 
DOUBLE  Pr<£ClSI3N  SCaLE(12),  SQRT2,  5J.  S»  START 
DlMENSlS.^  T T (224*  12) 

DIuENSld*  IT(224*12) 

DI  vtENSl  3n  SPNU6128) 

EQUIVALENCE  (IT(l).TT(D) 

Data  L»U»/224/.  lSCaLE/12/.  FNP/12.ODOO/ 

Pl/U.  141S927D00/.F0/20000.OD0U/.F3A/440.OD0O/ 

DEL  1/3. 00001 0000/. NP/1 2/ 

NBLnS/5/»LN6TH/288/»  MULTP/36/ 

RMS/1. 30448 10D00/.SQRT2/ 1.4 14^ 136000/. TP/O, 16128D00/ 
9NEJ3/1. 12201800C/. THROB/ 1 . 4 1 2337003/. TENDB/3, 1 622777D00/ 
NS«i>/95/ 

1 
2 


data 
data 
data 
data 
data 
data 
REwIno 
REWIND 


200 


100 

200 


C 

c 

c 


REWIND  3 
PAUSE 

IFjNSKP  »E0.  0)  60TB 
09  1 DO  I * 1*NSKP 
RcAD(I)  *»N 
OdiPHFsl. 122013«*1.5 
TW0D3=1. a22316»*2 
F30B=l.  U2D18*«3 
F40B=TW3U3**2 
F5qB=1.  U201S#»5 
F603  = TrlRD3*#2 
Fl2DB=F6«3**2 
F10DB=F6m3**3 
F24D3  = F103»#* 

F30DB=TEnd3»»3 
FBrMAT(5L2XI 18)/) 

change  thu  f3«  change  in  l@ut 

S/\  S PEAK  3|3NAL/SqRT(2.0)/RmS  N0ISE 
S/N  RAT  I 0^ 

SCaLE(1)=D.O 

SCaLE(2)=R^S*S0RT2/F6DB 

SCALE  (3)sR\*5*5QRT2/F4DB 

SCaLE(4)sRw5*SQRT2/F3DB 

SCaLE(5)*Rw5*SQRT2/TW3DB 

SCalE(6)srmS»SQRT2/0NEDB 

SCalE(7)*Rv!S»SQRT2 

SCalE(8)srwS*SORT2»0NEDB 


611 


SCALE (9)srmS*SQRT2*TN3DB 
SCALE  ( 10)  =Rv|S*S0RT2#F3DB 
SCALE  ( 1 1>=Rv1S*SQRT2#F4DB 
SCaLE(12)=RMS*SQRT2#F6DB 
MLTM1  sMJLTP-1 

LTt  sLSCALE»L0uT 

TW0  pi  =2. 0*P I 

FI  s3.5»FBW/Tp 

NUmBR  sMJLTP#LNGTH 

c change  this  fsr  change  in  l0ut 

NUmBR  = *5128 
FNuMBR  sNJMBR 

T.NUM3R  =(FNUMBR*1.0)/2.0 
0SCL  s 3*,3»F0*F0/(FNP«FB*) 

c oigItal  clock  frequency  2.424  megacycles 
TSCALE  sOELT*eSCL 

L2  =lTT 

L6  *LTT 


AA  vr*^  — * • T ] t 
LUivTiL/lI'y  i 14  XL. 


1 

1 


i 


112 


0®  60  Nsa,  V3LKS  /,  . 

reaou)  *on  cv..riL)t 

REWIND  2 
RR!TE<2)  SON 
08  55  Lsi.LSCALE 
IF(L.EQ.*>33  t*  23 
REWIND  2 
READ(2>  oosj 

C 19.2  KVDa  RANGE  (CHANGE  F8LL*WlNG  CARO  F8R  DIFF  RB  BIN) 

IF(M«0CN.3).NE.  1)  G8  T 8 30 
23  08  25  Jsa.NJNBR 
SJ  a J 

TX  a DELI.(SJ-TNUMBr) 

c geheraat*3n  sf  linear  fm  pulse 

S = DC8S(TW3PI»TX*(F0*FI*TX) ) 

25  SPn(J)  = SPN(J)*SCAlE(L)*S 
30  C9NTINUE 

03  55  I = A»  L3UT 

C C9VPUTE  STARl  = FIRST  Pe SITIVE  ZER8  CROSSING  rtITHiN  SAMPLE  INTERVAL 
Js72*I-7a 

40  IF(SPN(J>  >41. 42*43 

41  J = J ♦ 1 

I F ( SPN ( J A )41. 45.46 

42  J=J*1 

IF(SPN(Ji  >41,45.44 

43  Jaj+1 

G8  TS  40 

44  Jaj-i 

45  STARTaJ 
G9  T9  50 

46  P LlN  = SPN(J)/(SPN(J)-SPN(J-D) 

SJ  a j 

STaRT  = CUP  LIN-1. 5)#P  LIN  + 0.5)»U. 556-1. 0)*P  LIN  + SJ 
C C8UNT  nP  P8SiTlVE  ZER3  cR®SSlNGS  AND  THEN  CGMPUTE  ELAPSED  TIME  TT(I,L> 

50  D9  53  ICK3S3=l.NP 

51  J — J ♦ 1 

IF (SPN ( J)  )52.52»5i 

52  J-J*l 

IF(SPN(J>  >52.53.53 

53  C9nT!\UE 
SJ  a J 

TT(I.l>  s sj-start 
I F (SPN ( J>  >54.55.54 

54  P LlNaSPrt(j)/(SPN(J)-SPN(J-l)) 


64 


TT(I.L>  a(((P  L I N- a . 5 > *P  L I N*0. 5) »0. 566-; . 0) *P  LIN*TT(I.L) 
55  C9v)TlNUE 

M0pUL3  3*  C3UNTER  ©uTPuT  WITH  SLID»NG  PRESET 
D3  6"  JK-l.LSCALE 
ISET  *0 

0®  64  I*t,l9UT 

Change  thi©  fsr  change  in  l®ut 

IF  (MgD(lt7)*0  ,E0.  0)  ISET  = ISETM 
I T ( l » JK  > a M9D(IDINT(TSCALE*TT<I. J*> > ♦ 1 SET , 32) +0 
WRlTE(J)  IT 
WRlTE(6,t»l  i>  N 
60  C3MTINUE 


REWIND  1 
EN0  FILE 
REwInd  3 
ST®P 
END 


vwm  i I>ill  4 I 1/ 


!AL 
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P90GRAMS  L3 AOED 


C39E  LIMITS 


MAInI 

00600 

m 

01750 

E 1 61 

01751 

- 

47015 

47016 

- 

47037 

E 1 66 

47040 

- 

47043 

47044 

m 

4721  1 

DC9S 

47212 

m 

47231 

47232 

m 

47250 

DSIm 

47251 

m 

47254 

47255 

- 

47421 

47422 

- 

47472 

dlbg 

47473 

- 

47607 

DEXP 

47610 

«* 

47656 

47657 

- 

50026 

DM90 

50027 

• 

50072 

50073 

- 

50 1 22 

50123 

- 

50 1 24 

DA3s 

50125 

. 

50141 

A 1 62 

50142 

- 

50  1 56 

S 1 62 

50157 

• 

50i60 

50161 

- 

50200 

M I 62 

50201 

- 

50202 

50203 

- 

50217 

0162 

50220 

- 

50221 

50222 

m 

50242 

dble 

50243 

- 

50246 

50247 

- 

50260 

EI21 

50261 

m 

50277 

El  22 

50300 

- 

50301 

50302 

- 

50412 

AL9G 

50413 

• 

50425 

50426 

- 

50523 

EXP 

50524 

_ 

50537 

50540 

- 

50650 

M9D 

50651 

- 

50673 

50674 

- 

50720 

sngl 

50721 

50721 

50722 

- 

50770 

MC  | 6 1 

50771 

m 

51011 

MC|62 

51012 

m 

51040 

AS  I 22 

51041 

• 

51«  70 

M 1 22 

51171 

. 

51262 

0122 

51263 

- 

51361 

N * 22 

51362 

- 

5I366 

C I 2 1 
C I 7 1 
C I 8 i 

aint 

DINT 
I AT 
A I 7 1 
Cl  26 
Z I 80 
A I 8 i 
AS  I 66 
M I 66 
0166 
Nl  66 
C I 62 
C I 6 1 
Cl  16 
C 1 12 


51367 

51415 

51424 

51434 

51463 

51521 

51557 

51575 

51637 

51664 

51715 

52144 

52342 

52550 

52561 

52615 

52651 

52677 


5UM 
51423 
51433 
51462 
51520 
51556 
51574 
51636 
51663 
51714 
52143 
52341 
52547 
52560 
52614 
52650 
52676 
52743 
llL 


G'JriHOENTIAL 


VViU  j it\l 


1 


Z 1 To 

52744 

52761 

CIVT 

52762 

m 

53035 

ARU 

53036 

m 

53073 

CXlT/lsi 

53074 

- 

53075 

FtRTl/9 

53076 

- 

56616 

PIER 

56617 

• 

56703 

IPA 

56704 

- 

56721 

F9RT/MS-1 

56722 

- 

57357 

CONFIDENTIAL 


| RO 
1 2 1 
INE 
El  22 
C I 26 
E 1 2 1 
E 1 6 i 
Ml66 


53161 

55625 

55704 

50304 

51575 

50263 

47020 

52144 


i 


0166 

52342 

¥162 

50205 

C>  16 

52651 

A 166 

51736 

0162 

50224 

| NR 

53166 

M0D 

50676 

S 1 ®6 

51715 

oces 

47234 

A 162 

50144 

C 1 62 

52561 

S>  22 

51041 

0122 

51263 

SI62 

50163 

IOInT 

50737 

1 IX 

53652 

I NO 

53714 

IEN 

56531 

1ST 

53074 

1 AT 

51521 

E 1 66 

47046 

OLSg 

47475 

OEXp 

47661 

DM00 

50075 

FIEr 

566 1 7 

C 1 6 1 

52615 

OS  1 N 

47257 

DABS 

50127 

Nl  66 

52550 

C 1 8 1 

51424 

ZI80 

51637 

A 1 8 1 

51664 

Dl*T 

51463 

oble 

50251 

Cl  12 

52677 

ALSO 

50430 

¥122 

51171 

exp 

50542 

A I NT 

51434 

C 1 2 1 

51367 

C 1 7 1 

51415 

2 1 7 0 

52744 

A 1 22 

51047 

N 1 22 

51362 

A 1 7 1 

51557 

mC  1 62 

51012 

¥CI61 

50771 

gRRftR 

56631 

I 
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A*Ll 

ClVTARGll 
aROli 
ci  VTAQ 
ClVTQ 
UFA 
USN|  | 

|BA 
I IF 

no 

I IE 
I OE 
I EC 
I ef 

ISReT I 

1 ®tbl • 

WRITE l 
SETE6F 
rEadi 

B ACkS°  ACE » 
ENDFILE I 
rEwjnoi 


53036 

53020 

53061 

53000 

52762 

51533 

74506 

56524 

53646 

53641 

53634 

53103 

53076 

53142 

53341 

57212 

56740 

56722 

56735 

56727 

56732 

56743 


COiX'HDtNTIAL 
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co  co 
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o u 


o o 

o o 

uv  <o 

00  00 

o o 

o o 

O o 

u tp 


o o 

o o 

(vi  m 

<J»  Ov 

s § 

O o 

u u 


(M  <M 
-i  < 


(VJ  » 

o * 


1 

< 

h- 

Z 

4 

» 

< 

X 

_J 

< 

u 

\ 

I- 

r 

r> 

• 

ft 

3: 

1 

M 

z 

CM 

UJ 

(VI 

< 

z 

< 

< 

z 

z 

» 

II 

< 

n 

X 

w 

u 

tu 

X 

it 

(VJ 

3 

li- 

o 

a 

1 

UJ 

< 

< 

Z 

eu 

3 

r-t 

j 

z 

X 

X 
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